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Abstract We attach the degenerate signature (n,0, 1) to the projec- 
tivized dual Grassmann algebra P(/\ IR*^"+^)*) to obtain the Clifford algebra 
P(R* Q i) and explore its use as a model for euclidean geometry. We avoid 
problems with the degenerate metric by constructing an algebra isomorphism 
J between the Grassmann algebra and its dual that yields non-metric meet 
and join operators. We focus on the cases of n = 2 and n = 3 in detail, enu- 
merating the geometric products between k- and ^blades. We establish that 
sandwich operators of the form X — ^ gXg provide all euclidean isometries, 
both direct and indirect. We locate the spin group, a double cover of the di- 
rect euclidean group, inside the even subalgebra of the Clifford algebra, and 
provide a simple algorithm for calculating the logarithm of such elements. We 
conclude with an elementary account of euclidean rigid body motion within 
this framework. 



1 Introduction 



The work presented here was motivated by the desire to integrate the work 
of Study on dual quaternions ([Stu03]) into a Clifford algebra setting. The 
following exposition introduces the modern mathematical structures - pro- 
jective space, exterior algebra, and Cayley-Klein metrics - required to imbed 
the dual quaternions as the even subalgebra of a particular Clifford algebra, 
and shows how the result can be applied to euclidean geometry, kinematics, 
and dynamics. This paper is an extended version of a paper with the same 
name which was published in [Gunlla]. To help identify extensions, text not 
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contained in the original article appears dark blue. The original publication 
is available at www. springerlink . com. 



2 The Grassmann Algebra(s) of Projective Space 

Real projective n-space MP" is obtained from the (n + l)-dimensional 
euclidean vector space by introducing an equivalence relation on vectors 
X, y G \ {0} defined by: x ~ y x = Ay for some X ^ 0. That is, 

points in KP" correspond to lines through the origin in 

Grassmann algebra. The Grassmann, or exterior, algebra is 
generated by the outer (or exterior) product A applied to the vectors of M". 
The outer product is an alternating, bilinear operation. The product of a k- 
and m-vector is a (fc + m)-vector when the operands are linearly indepen- 
dent subspaces. An element that can be represented as a wedge product of 
k 1-vectors is called a simple fc- vector, or k-blade. The fc-blades generate the 
vector subspace /\'^(K"), whose elements are said to have grade k. This sub- 
space has dimension (^), hence the total dimension of the exterior algebra 
is 2". is one-dimensional, generated by a single element I sometimes 

called the pseudo-scalar. 

Simple and non-simple vectors. A fc-blade represents the subspace of 
K" spanned by the k vectors which define it. Hence, the exterior algebra 
contains within it a representation of the subspace lattice of M". For n > 3 
there are also fc-vectors which are not blades and do not represent a subspace 
of M". Such vectors occur as bivectors when n — A and play an important 
role in the discussion of kinematics and dynamics see Sect. 6. 

Projectivized exterior algebra. The exterior algebra can be projec- 
tivized using the same process defined above for the construction of MP" 
from K"'+-'^, but applied to the vector spaces /\*^(K""'"-'^). This yields the pro- 
jectivized exterior algebra W := P(A(K"+i)). The operations of 
carry over to W, since, roughly speaking: "Projectivization commutes with 
outer product." The difference lies in how the elements and operations are 
projectively interpreted. The fc-blades of W correspond to (fc — l)-dimensional 
subspaces of MP". All multiples of the same fc-blade represent the same pro- 
jective subspace, and differ only by intensity ([Whi98], §16-17). 1-blades cor- 
respond to points; 2-blades to lines; 3-blades to planes, etc. 

Dual exterior algebra. The dual algebra W* P(/\ ]R("+^'*) is formed 
by projectivizing the exterior algebra of the dual vector space (M""'""'^)* . De- 
tails can be found in the excellent Wikipedia article [Wik] , based on [Bou89] . 
W* is the alternating algebra of fc-multilinear forms, and is naturally iso- 
morphic to W; again, the difference lies in how the elements and operations 
are interpreted. Like W, W* represents the subspace structure of MP", but 
turned on its head: 1-blades represent projective hyperplanes, while n-blades 
represent projective points. The outer product aAb corresponds to the meet 
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rather than join operator. In order to distinguish the two outer products of 
W and W* , we write the outer product in W as V, and leave the outer prod- 
uct in W* as A. These symbols match closely the affiliated operations of join 
(union U) and meet (intersection fl), resp. 



2.1 Remarks on homogeneous coordinates 

We use the terms homogeneous model and projective model interchangeably, 
to denote the projectivized version of Grassmann (and, later, Clifford) alge- 
bra. 

The projective model allows a certain freedom in specifying results within 
the algebra. In particular, when the calculated quantity is a subspace, then 
the answer is only defined up to a non-zero scalar multiple. In some literature, 
this fact is represented by always surrounding an expression x in square 
brackets [x] when one means "the projective element corresponding to the 
vector space element x" . We do not adhere to this level of rigor here, since 
in most cases the intention is clear. 

Some of the formulas introduced below take on a simpler form which take 
advantage of this freedom, but they may appear unfamiliar to those used 
to working in the more strict vector-space environment. On the other hand, 
when the discussion later turns to kinematics and dynamics, then this pro- 
jective equivalence is no longer strictly valid. Different representatives of the 
same subspace represent weaker or stronger instances of a velocity or momen- 
tum (to mention two possibilities). In such situations terms such as weighted 
point or "point with intensity" will be used. See [Whi98], Book III, Ch. 4. 



2.2 Equal rights for W and W* 

From the point of view of representing MP", W and W* are equivalent. There 
is no a priori reason to prefer one to the other. Every geometric element in 
one algebra occurs in the other, and any configuration in one algebra has a 
dual configuration in the other obtained by applying the Principle of Duality 
[Cox87], to the configuration. We refer to as a point-based, and W* as a 
plane-based, algebra.^ 

Depending on the context, one or the other of the two algebras may be 
more useful. Here are some examples: 



^ We prefer the dimension-dependent formulation plane-based to the more precise 
hyperplane-based. We also prefer not to refer to the plane-based algebra as the dual 
algebra, since this formulation depends on the accident that the original algebra is 
interpreted as point-based. 



4 



Charles Gunn 



• Joins and meets. W is the natural choice to calculate subspace joins, 
and W* , to calculate subspace meets. See Sect. 2.2. 

• Spears and axes. Lines appear in two aspects: as spears (bivectors in 
W) and axes (bivectors in W*). See Sect. 2.2. 

• Euclidean geometry. W* is the correct choice to use for modeling eu- 
clidean geometry. See Sect. 3.2. 

• Reflections in planes. W* has advantages for kinematics, since it natu- 
rally allows building up rotations as products of reflections in planes. See 
Sect. 4.2. 

Bases and isomorphisms for W and W* . Our treatment differs from 
other approaches (for example, Grassmann-Cayley algebras) in explicitly 
maintaining both algebras on an equal footing rather than expressing the 
wedge product in one in terms of the wedge product of the other (as in 
the Grassman-Cayley shuffle product) ([Sel05], [Per09]). To switch back and 
forth between the two algebras, we construct an algebra isomorphism that, 
given an element of one algebra, produces the element of the second algebra 
which corresponds to the same geometric entity of EP". We show how this 
works for the case of interest n = 3. 

The isomorphism J. Each weighted subspace S of M.P^ corresponds to 
a unique element Sw of W and to a unique element Sw* of W* . We seek a 
bijection 3 : W ^ W* such that J{Sw) = Sw ■ If we have found J for the 
basis A:-blades, then it extends by linearity to multivectors. To that end, we 
introduce a basis for K'* and extend it to a basis for W and W* so that J 
takes a particularly simple form. Refer to Fig. 1. 

The canonical basis. A basis {e", e^, e^, e'^} of E'* corresponds to a 
coordinate tetrahedron for EP^, with corners occupied by the basis elements^. 
Use the same names to identify the elements of P(/\^(E'')) which correspond 
to these projective points. Further, let l" :— e° A e^ A e^ A e'^ be the basis 



Fig. 1 Fundamental 
tetrahedron with dual 
labeling. Entities in W 
have superscripts; entities 
in W* have subscripts. 
Planes are identified by 
labeled angles of two 
spanning lines. A rep- 
resentative sampling of 
equivalent elements is 
shown. 




^ We use superscripts for W and subscripts for W* since W* will be the more 
important algebra for our purposes. 
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element of P{/\'^{R'^}), and 1° be the basis element for P{/\°{R'^)). Let the 
basis for P(A (K^)) be given by the six edges of the tetrahedron: 

{e"i,e02,e03,ei2,e3i,e23} 

where e*-' :— e' A represents the oriented line joining e* and e-'.'^ Finally, 
choose a basis {E°,E\E2,E3} for P{/\^(R'^)) satisfying the condition that 
e* V E' =1°. This corresponds to choosing the z*'' basis 3- vector to be the 
plane opposite the basis 1-vector in the fundamental tetrahedron, oriented 
in a consistent way. 

We repeat the process for the algebra W* , writing indices as subscripts. 
Choose the basis 1-vector of W* to represent the same plane as E*. That 
is, J(E') = e^. Let Iq := Gq A ei A 62 A 63 be the pseudoscalar of the algebra. 
Construct bases for grade-0, grade-2, and grade-3 using the same rules as 
above for W (i. e., replacing subscripts by superscripts). The results are 
represented in Table 2.2. 



feature 


W 


W* 


0- vector 


scalar 1*^ 


scalar lo 


vector 


point {e*} 


plane {e^} 


bivector 


"spear" {e'^} 


"axis" {sij} 


trivector 


plane {E^} 


point {Ei} 


4- vector 




lo 


outer product 


join V 


meet A 



Table 1 Comparison of W and W* . 

Given this choice of bases for W and W* , examination of Fig. 1 makes 
clear that, on the basis elements, J takes the following simple form: 

3{e'):=E,, J{E') -.^ e,, 3{e'^) := e^i (1) 

where in the last equation, (ijkl) an even permutation of (0123). 

Furthermore, J(l°) = lo and J(I°) = lo since these grades are one- 
dimensional. To sum up: the map J is grade-reversing and, considered as 
a map of coordinate-tuples, it is the identity map on all grades except for 
bivectors. What happens for bivectors? In W, consider e°^, the joining line of 
points e° and (refer to Fig. 1). In W* , the same line is 623, the intersection 
of the only two planes which contain both of these points, 62 and 63. On a 
general bivector J takes the form: 

J(aoie"i -t- ao2e°2 -I- 0036°^ + ai2e^^ + 0316^^ + a23e^^) = 
023601 + 031602 + 012603 + 003612 + 002631 + 001623 

^ Note that the orientation of e^^ is reversed; this is traditional since Pliicker intro- 
duced these line coordinates. 
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The coordinate-tuple is reversed! See Fig. 2. Since J^^ is obtained from 
the definition of J by swapping superscripts and subscripts, we can consider 
3 : W ^ W* as a defined on both algebras, with the identity. The 
full significance of J will only become evident after metrics are introduced 
(Sect. 3.3). Appendix 1 contains a detailed description of how J is constructed 
in the n-dimensional case, and its close relation to the use of a nondegenerate 
metric to achieve the same goal. We now show how to use J to define meet 
and join operators valid for both W and W* . 

Projective join and meet. Knowledge of J allows equal access to join 
and meet operations. We define a meet operation A for two blades A,B € W: 



and extend by linearity to the whole algebra. There is a similar expression 
for the join V operation for two blades A,B G W*: 



We turn now to another feature highlighting the importance of maintaining 
W and W* as equal citizens. 

There are no lines, only spears and axes! Given two points x and y 
£ W , the condition that a third point z lies in the subspace spanned by the 
2-blade 1 := x V y is that x V y V z = 0, which implies that z = ax + f3y for 
some a, /3 not both zero. In projective geometry, such a set is called a point 
range. We prefer the more colorful term spear. Dually, given two planes x 
and y £ W* , the condition that a third plane z passes through the subspace 
spanned by the 2-blade 1 := xAy is that z = ax+/3y. In projective geometry, 
such a set is called a plane pencil. We prefer the more colorful term axis. 



AA B ^ J{J{A) AJ{B)) 



(2) 



AV B := J(J(A) V J(B)) 



(3) 



Fig. 2 A line in its dual 
nature as spear, or point 
range; and as axis, or 
plane pencil. 





Within the context of W and W*, lines exist only in one of these two 
aspects: of spear - as bivector in W - and axis - as bivector in W* . This 
naturally generalizes to non-simple bivectors: there are point-wise bivectors 
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(in W), and plane- wise bivectors (in W* .) Many of the important operators 
of geometry and dynamics wc will meet below, such as the polarity on the 
metric quadric (Sect. 3.1), and the inertia tensor of a rigid body (Sect. 6.2.2), 
map {W)2 to {W*}2 and hence map spears to axes and vice- versa. Having 
both algebras on hand preserves the qualitative difference between these dual 
aspects of the generic term "line" . 

We now proceed to describe how to introduce metric relations. 



3 Clifford algebra for euclidean geometry 

The outer product is anti-symmetric, so x A x = 0. However, in geometry 
there are important bilinear products which are symmetric. Wc introduce 
a real- valued inner product on pairs of vectors x • y which is a real-valued 
symmetric bilinear map. Then, the geometric product on 1-vectors is defined 
as the sum of the inner and outer products: 

xy := X • y -h X A y 

How this definition can be extended to the full exterior algebra is described 
elsewhere ([DFM09], [HS87]). The resulting algebraic structure is called a 
real Clifford algebra. It is fully determined by its signature, which describes 
the inner product structure. The signature is a triple of integers (p, n, z) 
where p + n + z is the dimension of the underlying vector space, and p, 
n, and z are the numbers of positive, negative, and zero entries along the 
diagonal of the quadratic form representing the inner product. We denote 
the corresponding Clifford algebra constructed on the point-based Grassmann 
algebra as P{Rp^n,z)', that based on the plane-based Grassmann algebra, as 

The discovery and application of signatures to create different sorts of 
metric spaces within projective space goes back to a technique invented by 
Arthur Caylcy and developed by Felix Klein [Kle49]. The so-called Cayley- 
Klein construction provides models of the three standard metric geometries 
(hyperbolic, elliptic, and euclidean) - along with many others! - within pro- 
jective space. This work provides the mathematical foundation for the inner 
prodiict as it appears within the homogeneous model of Clifford algebra. Since 
the Cayley-Klcin construction for euclidean space has some subtle points, is 
relatively sparsely represented in the current literature, and is crucial to what 
follows, we describe it below. 
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3.1 The Cayley-Klein Construction 

For simplicity we focus on the case n = 3. To obtain metric spaces inside 
M.P^ begin with a symmetric bihnear form Q on M.'^. The quadric surface 
associated to Q is then defined to be the points {x | (5(x, x) — 0}. For 
nondegenerate Q, a distance between points A and B can be defined by 
considering the cross ratio of the four points A, i?, and the two intersections 
of the hne AB with Q. Such a Q is characterized by its signature, there 
are two cases of interest for n — 3: (4, 0, 0) yielding elliptic geometry and 
(3,1,0) yielding hyperbolic geometry. These are point-based metrics; they 
induce an inner product on planes, which one can show is identical to the 
original signatures. By interpolating between these two cases, one is led to 
the degenerate case in which the quadric surface collapses to a plane, or to a 
point. In the first case, one obtains euclidean geometry; the plane is called the 
ideal plane. The signature breaks into two parts: for points, it's (1,0,3) and 
for planes it's (3, 0, 1). The distance function for euclidean geometry is based 
on a related limiting process. For details see [Kle49] or [Gunllb]. Warning: 
in the projective model, the signature (n, 0, 0) is called the elliptic metric, 
and euclidean metric refers to these degenerate signatures. 

Details of the limiting process. We restrict attention to a special family 
of Q parametrized by a real parameter e, and define an inner product between 
points X = (.To, Xi,X2, xs and Y = (yo, J/i, 2/2, ya) as follows: 



The inner product is not well-defined on M.P^ since different choices of 
representatives for the arguments will yield multiples of the inner product. 
But the expression 



is homogeneous in its arguments, hence well-defined (assuming neither argu- 
ment is a null vector). It is this normalized inner product which appears in 
the metric formulae below. 

gives the inner product for the signature (4,0,0) (elliptic space), and 
for (3,1,0) (hyperbolic space). For brevity we write these two inner 
products as (,)+ and (,)-, resp. 

The inner products above are defined on the points of space; there is an 
induced inner product on the planes of space formally defined as the adjoint 
operator of the operator Q^. One can easily show that for two planes u and 
V this takes the form : 



(X, Z)e := (3e(X, Y) := exoj/o + xiyi + X2J/2 + 2:32/3 



(4) 




(5) 



(X,X),(Y,Y), 



(u, v)^ := Q*{u,y) e ^Uot^o + "i^i +"2^2 + 

~ uqVo + e{uivi + U2V2 + U3V3) 



(6) 
(7) 
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where the ~ symbol denotes projective equivalence obtained by multiplying 
the first equation by e 7^ 0. For e £ {1,-1} the resulting signatures are 
equivalent to the point signatures; we can use the same notation for both 
points and planes in the formulae below. 

The metric quadric for elliptic space is {x | (x, x)+ = 0}. There are no real 
solutions; this is called a totally imaginary quadric. The points x | (x, x)+ > 
constitute the projective model of elliptic space: all of M.P^. The metric 
quadric for hyperbolic space is {x | (x,x)_ = 0}, the unit sphere in eu- 
clidean space. The points x | (x, x)_ < constitute the projective model of 
hyperbolic space: the interior of the euclidean unit ball. The sphere itself is 
sometimes called the ideal sphere of hyperbolic space. 

The euclidean metric is a limiting case of the above family of metrics as e 
grows larger and larger, from either the positive or negative side. For lim£_>oo, 
we arrive at the euclidean metric. Due to the homogeneity of (5) we allow 
ourselves to apply arbitrary non-zero scaling factors to the inner product. 
For points the limiting process yields: 

lim (X, Y)e = lim {exoyo + xiyi + X2V2 + xsys) 

~ lim {xoyo + - {xiyi + X2y2 + a^aya)) 
= xoVo 

where the second equation results from the first by multiplication by ^ . This 
is the signature (1, 0, 3). For planes u and v, one works with the adjoint form: 

lim ((u, v),) = lim (wqI^o + e(tiii'i + + "3^3)) 

e— ^C30 €— ^00 

~ lim {-unvo + uiwi + U2V2 + U3V3) 

= UlVl + U2V2 + U3W3 

which corresponds to the signature (3,0, 1). The point- and plane-signatures 
for euclidean space are complementary! 

The projective model of euclidean space consists of RP^ with the plane 
Xo — removed, the so-called ideal plane. For the points of this plane, the 
original inner product remains valid. These are euclidean free vectors, see 
discussion below. 

Distance and Angle Formulae. Let X and Y be two points, and Ixy 
their joining line. In general Ixy will have two (possibly imaginary) inter- 
section points Qi and Q2 with the metric quadric. The original definition 
of the distance of two points X and Y in these noneuclidean spaces relied 
on logarithm of the cross ratio of the points X, Y, Qi,, and Q2 ([Cox78]). 
By straightforward functional identities these formulae can be brought into 
alternative form. The distance d between two points X and Y in the elliptic 
(resp. hyperbolic) metric is then given by: 
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d = cosh-\— ^ ' ) (9) 

The familiar euclidean distance between two (non-homogeneous) points 

d = \J (xi - + (a-2 - y2f + (2:3 - ViY 

can be derived by parametrizing the above formulas by epsilon, and evalu- 
ating a carefully chosen limit as as e —t- 00. To simplify the treatment, we 
work with a projective line (n = 1), and introduce the family of inner prod- 
ucts parametrized by the real parameter e as above. We want to produce a 
euclidean metric on this line in which the basis vector ei is the ideal point. 
We consider two points X = {xo,xi) and Y = (yo,yi)- We can assume that 
a;o 7^ and j/o 7^ 0, since ei is the ideal point and doesn't belong to the 
euclidean line. Choose the projective representative so that xq — yo — 1- 
Then: 

(X, Y)e := exoyo + xiyi = e + xiyi 
By (8), the distance function associated to is determined by: 

cosd,(X,Y) = ^ (10) 

y((X,X),(Y,Y),) 

Abbreviate c?c(X, Y) as de- Now consider the limit as e — > 00. It's clear from 
(10) that lime_s.oo cosd^ = 1, that is, the distance goes to zero. Define a new 
distance function de '■= \fld^. The scaling factor -^/e prevents the distance 
from going to zero in the limit. Instead one obtains: 

lim dl = (xi - y^Y (11) 

€—^00 

which is clearly equivalent to the euclidean distance between the two points 
X and Y. For the details consult Appendix 2. 

In all three geometries the angle a between two oriented planes u and v 
is given by (where (,) represents the appropriate inner product): 

a^cos-\-^^ ) 

'v/((u,u)(v,v))^ 

Polarity on the metric quadric. For a Q and a point P, define the 
set P-L := {X I (3(X,P) = 0}. When P-^ is a plane, it's called the jiolar 
plane of the point. For a plane a, there is also an associated polar point 
defined analogously using the "plane-based" metric. Points and planes with 
such polar partners are called regular. In the euclidean case, the polar plane 
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of every finite point is the ideal plane; the polar point of a finite plane is the 
ideal point in the normal direction to the plane. Ideal points and the ideal 
plane are not regular and have no polar partner. The polar plane of a point 
is important since it can be identified with the tangent space of the point 
when the metric space is considered as a differential manifold. Many of the 
peculiarities of euclidean geometry may be elegantly explained due to the 
degenerate form of the polarity operator. In the Clifford algebra setting, this 
polarity is implemented by multiplication by the pseudoscalar. 

Free vectors and the euclidean metric. As mentioned above, the tan- 
gent space at a point is the polar plane of the point. Every euclidean point 
shares the same polar plane, the ideal plane. In fact, the ideal points (points 
of the ideal plane) can be identified with euclidean free vectors. A model 
for euclidean geometry should handle both euclidean points and euclidean 
free vectors. This is complicated by the fact that free vectors have a natu- 
ral signature (3,0,0). However, since the limiting process (in Cayley-Klein) 
that led to the degenerate point metric (1,0,3) only effects the non-ideal 
points, it turns out that the original non-degenerate metric, restricted to the 
ideal plane, yields the desired signature (3, 0, 0). As we'll see in Sect. 4.1 and 
Sect. 5.2, the model presented here is capable of mirroring this subtle fact. 

To learn more about the mathematics, see [Kle49], [Cox87] [Cox78], and 
[Kow09]. For a modern formulation with a good collection of formulae see 
[TL97]. 



3.2 A model for euclidean geometry 

As noted above, the euclidean inner product has signature (1, 0, 3) on points 
and (3, 0, 1) on planes. If we attach the first signature to W, we have the 
following relations for the basis 1-vectors: 

(e")2 = 1; {e'f = [e'f = {e^f = 

It's easy to see that these relations imply that, for all basis trivectors £'i, 
Ef = 0. But the trivectors represent planes, and the signature for the plane- 
wise euclidean metric is (3,0,1), not (0,0,4). Hence, we cannot use W to 
arrive at euclidean space. If instead, we begin with W* , and attach the plane- 
wise signature (3, 0, 1), we obtain: 

(eo)2 = 0; (eO' = {e^f = {e^f = 1 

It is easy to check that this inner product, when extended to the higher grades, 
produces the proper behavior on the trivectors, since only Eq = eie2e3 has 
non-zero square, producing the point-wise signature (0,1,3) (equivalent to 
the signature (1,0,3)). Hence, W* is the correct choice for constructing a 
model of euclidean geometry. 
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Counterspace. What space does one obtain by attaching the signature 
(3,0, 1) to W? One obtains a different metric space, sometimes called polar- 
euclidean space or counterspace. Its metric quadric is a point along with all 
the planes passing through it (dual to the euclidean ideal plane and all the 
points lying in it)^. Like euclidean space, it arises as a limiting case of the 
Cayley-Klein construction, when one lets e — )■ 0. Instead of flattening out the 
metric quadric into a plane, this limit contracts it to a point. See [Con08], 
pp. 71ff., for a related discussion. 

We retain W, the point-based algebra, solely as a Grassmann algebra, 
primarily for calculating the join operator. All euclidean metric operations 
are carried out in W* . Or equivalently, we attach the metric (0,0,4) to W, 
forcing all inner products to zero. Due to the more prominent role of W* , 
the basis element for scalar and pseudoscalar in W* will be written without 
index as 1 and I; we may even omit 1 when writing scalars, as is common in 
the literature. 



3.3 3, metric polarity, and the regressive product 

We can now appreciate better the significance of J : T4^ — > W* . Consider the 
map n : W ^ W defined analogously to J in (1): 

n(e*) := E\ U{E') := e\ U{e'^) := e''^ (12) 

n is the same as J, but interpreted as a map to W instead of W* . It's easy 
to see that 11 is the polarity on the elliptic metric quadric with signature 
(4,0,0). Many authors (see [HS87]) define the meet operation between two 
blades A, B G W (also known since Grassmann as the regressive product) via 
n(n(>l) A 11(B)), where A is the exterior product in W. One can define a 
similar join operator in W* . We prefer to use J for this purpose (see (2)) since 
it provides a projective solution for a projective (incidence) problem, and it 
is useful on its own (see for example Sect. 6.2.2), while II, being a foreign 
entity, must always appear in the second power so that it has no side-effects. 
To distinguish the two approaches, we suggest calling II the metric polarity 
and J, the duality operator, consistent with mathematical literature. For an 
n-dimensional discussion and proof, see Appendix 1. 



* Blurring the distinction between these two spaces may have led some authors to 
incorrect conclusions about the homogeneous model, see [Li08], p. 11 
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Fig. 3 The basis 1- and 
2-vectors for P(K^,q_i). 
eo is the ideal line, Eg is 
the origin. 

4 The euclidean plane via P(]R* ^ ^) 

Due to the combination of unfamiliar concepts involved in the algebras 
P(M* oi) ~ notably the dual construction and the degenerate metric - we 
begin our study with the Clifford algebra for the euclidean plane: P(M2 o i)- 
Then, when we turn to the 3D case, we can focus on the special challenges 
which it presents, notably the existence of non-simple bivectors. A basis for 
the full algebra of P(IR2 o i) given by 

{1 := lo,eo,ei,e2,Eo := eie2,Ei := e2eo,E2 := eoei,I := 606162} 

with the relations {gq = 0; = 63 = 1}. See Fig. 3. 

Consequences of degeneracy. The pseudoscalar I satisfies — 0. 
Hence, is not defined. Many standard formulas of geometric algebra are, 
however, typically stated using ( [DFM09, HS87]), since that can simplify 
things for nondegenerate metrics. As explained in Sect. 2.1, many formulas 
remain projectively valid when I^^ is replaced by I; in such cases this is the 
solution we adopt. 

Notation. We denote 1-vectors with bold small letters, and 2-vectors with 
bold capital letters. We will use the term ideal to refer to geometric elements 
contained in projective space but not in euclidean space. Then eg is the ideal 
line of the plane, ei is the line x = and e2, the line y = 0. Eq is the origin 
(1,0,0) while El and E2 are the ideal points in the x— and y— direction, 
resp. Points and lines which are not ideal, are called finite, or euclidean. 

The natural embedding of a euclidean position x = (a;, y) we write as 
i(x) = Eq + .xEi -|-yE2. A euclidean vector v = (x, y) corresponds to an ideal 
point (see Sect. 3.1); we denote its embedding with the same symbol i(v) = 
xEi-|-yE2. We sometimes refer to such an element as a, free vector. Conversely, 
a bivector wEq + a;Ei + J/E2 with w ^ Q corresponds to the euclidean point 
( — , — ). We refer to w as the intensity or weight of the bivector. And, we 
write A to refer to i~^(A). The line ax + hy + c ~ maps to the 1-vector 
ceo + aei + he2- A line is euclidean if and only if + 6^ 7^ 0. 
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Table 2 Geometric product in P(]RJ q 1) 



The multiplication table is shown in Table 4. Inspection of the table reveals 
that the geometric product of a k- and ^-vector yields a product that involves 
at most two grades. When these two grades arc \k — l\ and k + l,we can write 
the geometric product for 2 arbitrary blades A and B as 

AB = A B + A A B 

where • is the generalized inner product, defined to be (AB)|fc_;| ([HS87]). The 

only exception is (l, k) = (2, 2) where the grades \k — l\ =0 and \k — l\+2 = 2 
occur. Following [HS87] we write the grade-2 part as: 

A X B := (AB)2 = ^(AB - BA) 

where A and B are bivectors. This is called the commutator product. Since 
all vectors in the algebra are blades, the above decompositions are valid for 
the product of any two vectors in our algebra. 



4-1 Enumeration of various products 

We want to spend a bit of time now investigating the various forms which the 
geometric product takes in this algebra. For this purpose, define two arbitrary 
1-vectors a and b and two arbitrary bivectors P and Q with 

a = aoco + oiei + 0262, etc. 

These coordinates arc of course not instrinsic but they can be useful in un- 
derstanding how the euclidcan metric is working in the various products. See 
the companion diagram in Fig. 4. 

1. Norms. It is often useful to normalize vectors to have a particular inten- 
sity. There are different definitions for each grade: 
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• 1-vectors. = a-a = ai+a|. Define the norTn of a to be II a| I := y/a. ■ a. 

Then -rr—rr is a vector with norm 1, defined for all vectors except eo and 

l|a|| 

its multiples. In particular, all euclidean lines can be normalized to have 
norm 1. Note that when a is normalized, then so is —a. These two lines 
represents opposite orientations of the line^. 

• 2-vectors. = P • P = PqEq = — pg- Define the norm of P to be po 

and write it ||P||. Note that this can take positive or negative values, 

p 

in contrast to \/P P Then TT— T is a bivector with norm 1 , defined for 

all bivectors except where po = 0, that is, ideal points. In particular, all 
euclidean points can be normalized to have norm 1. This is also known 
as dehomogenizing. 

• 3-vectors'. Define S : P{/\^R^*) P(A°M^*) by S(al) = al. This 
gives the scalar magnitude of a pseudoscalar in relation to the basis 
pseudoscalar I. Wc sometimes write ^(al) for the same. In a non- 
degenerate metric, the same can be achieved by multiplication by 

a — P 

2. Inverses, a"-'^ = and P""'^ = — . for euclidean a and P. 

a a P P 

3. Euclidean distance. For normalized P and Q, ||P V Q|| is the euclidean 

distance between P and Q. 

4. Free vectors. For an ideal point V (that is, a free vector) and any nor- 
maHzed euclidean point P, ||V||oo := ||VVP|| = ^yv1 + V2 is the length 

V 

of V. Then — — - — is normalized to have length 1. 

II V lloo 

5. a A P = {aoPo + o,iP2 + 02^2)1 vanishes only if a and P are incident. 
Otherwise, when a and P are normalized, it is equal to the signed distance 
of the point to the line times the pseudoscalar I . 

6. P • a = {p2ai — pia2)eo + poct2ei — po«ie2 is a line which passes through 
P and is perpendicular to a. Reversing the order changes the orientation 
of the line. 

7. a A b =: T is the intersection point of the lines a and b. For normalized a 
and b, ||Tj| = sina where a is the angle between the lines. Reversing the 
order reverses the orientation of the resulting point. 

8. a • b = cos a for normalized vectors a and b. Which of the two possible 
angles is being measured here depends on the orientation of the lines. 

9. P V Q is the joining line of P and Q.. 

10. P X Q =: T is the ideal point in the direction perpendicular to the direction 
of the line P V Q. 

11. al = fliEi + a2E2 is the polar point of the line a: the ideal point in the 
perpendicular direction to the line a. All lines parallel to a have the same 
polar point. 

12. PI — paeo is the polar line of the point P: for finite points, the ideal line, 
weighted by the intensity of P. Ideal points have no polar line. 



^ Orientation is an interesting topic which lies outside the scope of this article. 
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Fig. 4 A selection of the 
geometric products be- 
tween various fe-blades. 
Points and lines are as- 
sumed to be normalized. 
Ideal points are drawn as 
vectors, distances indi- 
cated by norms. 




13. = 0. This is equivalent to the degeneracy of the metric. Notice that this 
fact has no effect on the vaUdity of the above calculations.® 



4.1.1 Exercises 

Rather than presenting finished results based on the above formulas, we 
present the material in the form of exercises for the reader to work through. 
Readers who choose not to work through the results are still recommended 
to read through them, since what follows will build on the results of the 
exercises. Exercises marked with an asterix are more challenging. 

1. Angles between vectors. For normalized ideal points U and V, 

(WP) • (UVP) =cosa 

where a is the angle between the vectors. 

2. Projection onto a line. For a euclidean line a and euclidean point P, 
show that (P • a)a^^ is the orthogonal projection of the point P onto a. 

• What does (a • P)P~^ represent? 

• Compare (P • a)a and (P • a)a^^. 

3. For normalized arguments, show that ||a- P|| = 1. 

4. Parallel lines. What is the condition that a and b are parallel? In the 
expression for the distance of two lines above, what result is obtained when 
a and b are parallel? 

5. What changes have to be made to the above formulas when the arguments 
are not normalized? 



^ In fact, the validity of most of the above calculations requires that I2 = 0. 
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Fig. 5 The three triangle 
centers M (centroid), 
T (orthocenter), and P 
(circumcenter) lying on 
the Euler line e of the 
triangle ABC. 

6. Ideal elements. What changes have to be made to the above formulas 
when the arguments are not finite? In particular, are the formulas (6) and 
(7) above valid when one or both of P and Q are ideal points? 

7. P X Q = (PVQ)I. 

8. Show that for normalized A,B, and C, the area of AAHC is given by 
(AVB) AC. 

9. Triangle centers. Given a euclidean triangle ABC, this exercise shows 
how to use P(M2oi) to calculate the four classical triangle centers, and 
prove three lie on the Euler line. Consult Fig. 5. 

a. Show Z\ABC can be assumed to have normalized corners A = Eg, 
B = Eq + &iEi, and C = Eq + ciEi + C2E2. 

b. Define the edge-lines a := B V C, b := C V A, and c A V B. 

c. Show that the median niA, perpendicular bisector p^, altitude t^ and 
angle bisector n^ associated to the pair a and A are given by the fol- 
lowing formulas: 

• nia = (B + C) V A 

• Pa = (B + C)-a 

• ta = A • a 

• n^ = ||c||b + ||b||c 

d. Derive analogous formulas for the lines associated to b and c. 

e. Show that hIq, mf,, and nic are co-punctual; their common point is 
the centroid M of the triangle. [Hint: to show that the three lines go 
through the same point you can show the outer product of all three is 
0. If you work with the fully general forms for A, B, and C (don't use 
9a), you can also show that the expression for the intersection of two 
of the lines is symmetric in {a,b,c} and {A, B,C}.] 

f. Do the same for the triple of lines 

• (pa, pfc, Pc), to obtain the circumcenter P, 

• (ta, th, tc), to obtain the orthocenter T, and 

• (n^, n^, lie), to obtain the incenter N. 
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g. Show that M, P, and T he on a hne e (the Euler line of the triangle). 
[Hint: to show three points he on a hne, show that the join (V) of two 
of the points has vanishing outer product with the third. 

h. Show that M hes between P and T on the Euler line, and is twice as 
far from T as from P. 



4.-2 Euclidean isometries via sandwich operations 

One of the most powerful aspects of Clifford algebras for metric geometry is 
the ability to realize isometries as sandwich operations of the form 

X ^ gXg-i 

where X is any geometric element of the algebra and g is a specific geometric 
element, unique to the isometry. g is in general a versor, that is, it can be 
written as the product of 1-vectors ([HS87]. Let's explore whether this works 
in P(M^,o,i)- 

Reflections. Let a := ep — ei (the line x = 1), and P a normalized point 
Eq + xEi + yEi2- Simple geometric reasoning shows that refiection in the 
line a sends the point (x, y) to the point (2 — x, y). Let's evaluate the versor 
operator: 

P' := aPa^ = aPa = . . . 

= (eo - ei)(Eo + xEi + yE2)(eo - ei) 

= (eo - ei)(eie2 + a:e2eo + yeoei)(eo - ei) 

= (eo - ei)(eie2eo - eie2ei - a;e2eoei - yeoeiei) 

= (eo - ei)((l - a::)eoeie2 + e2 - ye^) 

= (eoe2 - (1 - a:)eieoeie2 - eie2 + j/eieo) 

= {x- 2)el - eO - ye2 

= eO + (2 - a;)el + ye2 

where the final step is obtained by dehomogenizing.^ 

This algebra element corresponds to the euclidean point {2 — x,y), so the 
sandwich operation is the desired refiection in the line a. We leave it as 
an exercise for the interested reader to carry out the same calculation for a 
general line. 

Direct Isometries. By well-known results in plane geometry. 

The result of carrying out reflections in two lines one after the other (the 
composition of the reflections) is a rotation around the intersection point of 



Without homogenizing, the orientation of the point is reversed, as it probably should 
be. 
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the two lines, through an angle equal to twice the angle between the two lines, 
unless the two lines are parallel, in which case the composition is a translation 
in the direction perpendicular to the two lines, through a distance equal to 
twice the distance between the two lines. 

Translating this into the language of the Clifford algebra, the composition 
of reflections in lines a and b will look like: 

P' = b(aPa)b 
= TFT 

where we write T ba, and T is the reversal of T. 

Note that the intersection of the two lines will be fixed by the resulting 
isometry. There are two cases: the point is ideal, or it is euclidean. In the 
case of an ideal point, the two lines are parallel and the composition is a 
translation. Let's look at an example. 

Retaining a as above, define the normalized line b :— 2eo — ei, the line 
X — 2. By simple geometric reasoning, the composition "reflect first in a, then 
in b" should be the translation {x,y) — > (a; + 2,y). Defining T := ba, the 
sandwich operator looks like: TFT. Calculate the product T = 1 — E2, and 
TFT = Eo + (2 + a;)Ei + E2. This shows that TFT is the desired translation 
operator. One can generalize the above to show that a translation by the 
vector (xo^yo) is given by the sandwich operation TFT where 

T:=l + i(j/oEi-xoE2) (13) 

It's interesting to note that TF and FT are both translations of {x, y) 
(a; + 1, y), so one doesn't need a sandwich to implement translations, but for 
simplicity of representation we continue to do so. 

Rotations. Similar remarks apply to rotations. A rotation around a nor- 
malized point R by an angle 9 is given by 

T = cos(^) + sin(^)R 

This can be checked by substituting into (13) and multiplying out. We'll ex- 
plore a method for constructing such rotators using the exponential function, 
in the next section. See [Gunllb] for a more detailed discussion including 
constructions of glide refiections and point reflections. 

4.2.1 Reflections in points and in hyperplanes 

When a is a 1-vector, then aXa is the reflection in a. In the dual algebra, a 
represents a hyperplane (in this case, a line); in the standard algebra, a point. 
Because reflections in planes are often more practically useful than reflections 
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in points, this can be seen as an advantage for the dual approach used here. 
In a non-degenerate metric this is not a significant advantage, because in such 
metrics a reflection in a hyperplane is also a reflection in the polar point of 
the plane; and vice-versa. 

This is also good place to point out that a common form for a reflection 
in a hyperplane, in the vector space model, includes a minus sign ([DFM09], 
p. 168). Thus, the reflection in the plane whose normal vector is P is given 
by: 

X ^ PXP 

This minus sign is due to the fact that the algebra is point-based but the 
desired reflection is in a hyperplane. Without minus sign, the expression PXP 
represents a reflection in the vector P. In 3 dimensions, this is a rotation 
of 180 degrees around the line AP. To obtain the reflection in the plane 
p orthogonal to P, one must compose the vector reflection with the point 
reflection in the origin, which is achieved by multiplying by —1. This undoes 
the rotation around the line AP and introduces a reflection in the orthogonal 
plane. This yields the expression —PXP as the form for a reflection in the 
plane orthogonal to P. Compare this to the homogeneous model presented 
here, where reflections in points and in planes are represented without any 
extra minus signs. 

4.2.2 Exercises 

1. A glide reflection is the product of a reflection in a line with a transla- 
tion parallel to the line. Describe how to realize a glide reflection using a 
sandwich operation. 

2. A point reflection in a point P is an isometry that sends each point Q to its 
reflected image on the "other side" of P. In particular, the image of Q lies 
on the line P V Q on the other side of P an equal distance to P. Show that 
Q —> PQP realizes this point reflection. How is this consistent with the 
theory developed above that such rotors are rotations? [Hint: cos ^ = 0.] 

3. *Show that the above formulas for the action of reflections, translations, 
and rotations on points are also valid when P is replaced by a euclidean 
line a. 



4-3 Spin group, Exponentials and Logarithms 

We have seen above in (13) that euclidean rotations and translations can be 
represented by sandwich operations in P(M.2 q i), in fact, in the even subal- 
gebra P(M;+ J. 
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Definition 1. The spin group Spin(2,0,l) consists of elements g of the 
even subalgebra P(IR21) i) such that gg = 1. An element of the spin group is 
called a rotor. Some authors ([Per09]) refer to an element of the spin group 
as a spinor but this conflicts with the accepted definition of spinor in the 
mathematics community, so we avoid using the term here. 

Write g = si + M where M = (g)2 = moEo + miEi + m2E2. Then gg = 
+ jtiq — 1. There are two cases. 

• mo 7^ 0, so M is a euclidean point. Then there exists 6* 7^ such that s = 

M 

cos (0) and niQ — sin (9), yielding g = cos {6) + sin (6')N where N 



sma 

is a normalized point, hence ~ —1. Thus, the formal exponential e*'^ 
can be evaluated to yield: 



e 



i=Q 



= cos {t) + sin (i)N 

Hence, the rotor g can be written as an exponential: g = e^'^. 
mo = 0, so M is an ideal point. Then we can assume that s = 1 (if 
s = — 1, take the element — g with the same sandwich behavior as g). 
Also, M = miEi + m2^2- Again, the formal exponential e*'^ can be 
evaluated to yield: 



= l + tM 

So in this case, too, g = has an exponential form. 
The above motivates the following definitions: 

Definition 2. A rotator is a rotor whose bivector part is a euclidean point. 
A translator is a rotor whose bivector part is a ideal point. 

Definition 3. The logarithm of a translator g = si + M e P(lR2,o 1) i*^ 

M 

smce e = g. 

Definition 4. Given a rotator g = sl+moEo+miEi+m2E2 € Spin(2,0, 1). 
Define 6 := tan^^(TOo:s) and N := j^^^- Then the logarithm of g is 6'N, 
smce e = g. 
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4.3.1 Exercises 

1. Show that a translation of the point P by the vector {xq, j/q) can be written 
as e'^Pe-^ for M = 1 + ^(yoEi - xqEz). 

2. Show that a rotation of the point P around the point {xo^yo) by angle 6 
can be written as e'^Pe"^ for M = sin (|)(Eo + xqEi + z/oE2). 

3. Deduce a condition on two rotors gi and g2 so that their product is a 
translator. 

Lie groups and Lie algebras. The above remarks provide a realization 
of the 2-diniensional cuclidean direct isometry group SE(2) and its Lie al- 
gebra se(2) within P(M2t) i)- The Spin group Spin(2,0, 1) forms a double 
cover of SE(2) since the rotors g and —g represent the same isometry. Within 
P(M2'o i): the spin group consists of elements of unit norm; the Lie algebra 
consists of the pure bivectors plus the zero element. The exponential map 
X — > e-^ maps the latter bijectively onto the former. This structure is com- 
pletely analogous to the way the unit quaternions sit inside P(IR3^g) and 
form a double cover of 50(3). The full group including indirect isometries is 
also naturally represented in P(K2 o i) ^^ic group generated by reflections 
in lines, sometimes called the Pin group. 



4-4 Guide to the literature 

There is a substantial literature on the four-dimensional even subalgebra 
P(K2'o i) with basis {1, Eq, Ei, E2}. In an ungraded setting, this structure 
is known as the planar quaternions. The original work appears to have been 
done by Study ([Stu91], [Stu03]); this was subsequently expanded and refined 
by Blaschke ([Bla38]). Study's parametrization of the full planar euclidean 
group as "quasi-elliptic" space is worthy of more attention. Modern accounts 
include [McC90]. ^ 



5 P(]R3 Q j^) and euclidean space 

The extension of the results in the previous section to the three-dimensional 
case P(IR3 p 1) is mostly straightforward. Many of the results can be carried 
over virtually unchanged. The main challenge is due to the existence of non- 
simple bivectors; in fact, most bivectors are not simple! (See Sect. 5.1 below.) 
This means that the geometric interpretation of a bivector is usually not a 
simple geometric entity, such as a spear or an axis, but a more general object 



* which however confuses the euclidean inner product on vectors with the inner 
product on points. 
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known in the classical literature as a linear line complex, or null system. Such 
entities are crucial in kinematics and dynamics; we'll discuss them below in 
more detail. 

Notation. As a basis for the full algebra we adopt the terminology for 

the exterior algebra W* in Sect. 2.2, interpreted as a plane-based algebra. We 
add an additional basis 1- vector satisfying = 1. eo now represents the ideal 
plane of space, the other basis vectors represent the coordinate planes. Eq is 
the origin of space while Ei = 606362 is the ideal point in the x-dircction, 
similarly for E2 and E3. The bivector eoi is the ideal line in the x = plane, 
and similarly for 602 and 603. 623,631, and 612 are the x-, y-, and 2;-axis, 
resp. We use i again to denote the embedding of euclidean points, lines, and 
planes, from RP^ into the Clifford algebra. 

We continue to denote 1-vectors with bold small Roman letters a; trivec- 
tors will be denoted with bold capital Roman letters P; and bivectors will be 
represented with bold capital Greek letters H.^ 

We leave the construction of a multiplication table as an exercise. Once 
again, most of the geometric products of two vectors obey the pattern AB = 
A ■ B + A A B. Two new exceptions involve the product of a bivector with 
another bivector, and with a trivector: 

H* = H • * + H X * + H A * (14) 
HP = S-P + HxP (15) 

Here, as before, the commutator product Ax B := ^{AB — BA). 

We now describe in more detail the nature of bivectors. We work in W* , 
since that is the foundation of the metric. As a result, even readers familiar 
with bivectors from a point-based perspective will probably benefit from going 
through the following plane-based development. 



5.1 Properties of Bivectors. 

We begin with a simple bivector H := a A b where a and b are two planes 
with coefficients {aj} and {bi}. The resulting bivector has coefficients 

Pij := aibj - ajbi {ij e {01, 02, 03, 12, 31, 23}) (16) 

These are the plane-based Pliicker coordinates for the intersection line (axis) 
of a and b. Clearly H A H = 0. Conversely, if 

S A S = 2(poiP23 + P02P31 + P03Pl2)I = 
for a bivector H, the bivector is simple [Hit03]. 



A convention apparently introduced by Klein, see [Kle72]. 
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Given a second axis = c A d, the condition H A 4> = implies they 
have a plane in common, or, equivalently, they have a point in common. For 
general bivectors, 

H A * = (poi'723 + 7502931 + P03912 + Puqos + Pziqm + P23goi)I (17) 

The parenthesized expression is called the Pliicker inner product of the two 
lines, and is written (H, #)p. With this inner product, the space of bivectors 
P{/\^{W^)*) is the Cayley-Klein space *B := P(M3,3)^ ^nd the space of hues 
is the quadric surface L"^'"^ = {H | (H, H)p = 0} C 03. *B is sometimes 
called Klein quadric. When (17) vanishes, the two bivectors are said to be in 
involution. 

Pencils of line complexes. Two bivectors H and $ span a line in *B, 
called a line complex pencil, or bivector pencil. Points on this line are of the 
form aS + for a, /3 € M, both not 0. Finding the simple bivectors on this 
line involves solving the following equation in the homogeneous coordinate 



This equation can be undetermined, or have 0, 1, or 2 real homogeneous 
roots. See [PWOl] for details. Finding the intersection of a line in 58 with 
L'^'^ is a common procedure in line geometry. It's used below to calculate the 
axis of a non-simple bivector. 

Null polarity associated to a non-simple bivector. For a fixed H, 
the orthogonal complement H-'- := | (H,^')p = 0} is a 4-dimensional 
hyperplane of *8 consisting of all bivectors in involution to H. The intersection 
with L^''^ is a 3-dimensional quadratic submanifold called a linear line 
complex. When H is simple, this consists of all lines which intersect H. 

Assume that H is not simple. Then H determines a coUineation, the har- 
monic homology Ha : MP^ — MP^, with center H and axis H-'" given by the 
following formula: 



Here (, ) is the canonical pairing of a vector and a dual vector in the vector 
space underlying MP^. When expressed in terms of the Pliicker inner 
product, this takes the second form given above. 

He is a collineation of 58 which fixes H and H-*- pointwise. It's a kind 
of mirror operation when applied to a point Let X be the line in 58 
joining H and and Hj^ the intersection of X with H-*-. Then H3($) is 
the harmonic partner of $ with respect to the fixed points H and Hj^, the 
point on X satisfying the cross ratio condition CR(0, H, H_l) — —1. If 



A = a : /3: 



= (aH + /3*, aH + /3*)p 
= a2(H, H)p + 2a/3(*, H)p + *)p 



H(*) := (H,H^)*-2(*,H^)H 
:= (H,H)p*-2(*,H)pH 



(18) 
(19) 
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$ = aS + /3Hj_, then = aS, — /3Hj_ (Exercise). When $ is simple, then 
so is (Exercise). To proceed, we need a theorem relating projectivities of 
M.P^ and projectivities of *B preserving L'^'^. See [PWOl], p. 144, for a proof. 

Theorem 1. The group of projectivities o/MP^ and the group of coUineations 
of *B preserving L'^'^ are isomorphic. 

This result confirms that the harmonic homology H3 has an induced ac- 
tion on MP^, the null polarity associated to H. An simple element of H-*-, 
considered as a line in M.P^, is called a null line of the null polarity. First of 
all, how do the null lines of H appear? Consider a point P, and define the null 
plane of P to be Ns(P) := H V P. Then iVH(P) V P 0, so P hes in its nuh 
plane. Also, let Q be another point of Ns{P)- Then apply the associativity 
of the V product to obtain: 



This shows that the line P V Q is a null line of H. Conversely, any null line of 
H passing through P can be written in the form P V Q, and one can reverse 
the reasoning to conclude it must lie in Ns{P)- 

One can also dualize the discussion to define the null point of a plane 
Ns{sl) := HA a. Ns,{Ns{sl)) = a by (36) , showing that Ns is an involution, 
hence in fact deserves the name null polarity. 

Comparison to versors. Sect. 5.5 shows that euclidean isometrics can 
be represented in 3D also with versor operators of the form gXg. These 
isometrics are projective transformations of M.P^ that preserve the metric 
quadric. The null polarity Ns is also a projective transformation of M.P^ but 
it is realized within the Clifford algebra in a different way^'^. The action on 
planes, lines, and points is given by: 



Of course, once the action on points, on lines, or on planes is given, the 
action on the other two types of elements is determined. One can construct 
other projective transformations by composing such null polarities. For exam- 
ple, if H and # are in involution, then the composition N^oNs = NsoN^ is 
an involution which preserves each null line common to both line complexes 
$ and H. By composing three such null polarities, one arrives at the polarity 
with respect to a regulus (the common null lines of all 3 line complexes) . For 
details see [Wei35]. 



= (iVH(P) V Q = (H V P) V Q 
= H V (P V Q) 



iVH(a) 
A^h(*) 
^s(P) 



S(HAH)*-2S(*AH)H 
EVP 



HA a 



(20) 
(21) 
(22) 



^'^ Actually it's more accurate to say it is realized in the Grassman algebra since it 
doesn't involve the inner product 
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By taking all possible compositions of these null polarities, one arrives at 
a representation of the full projective group of M.P^ within the Grassman 
algebra. These representations do not have the simplicity and homogeneity 
of the versor representation of isometries. Still, this representation of the 
projective group seems worth investigating. For example, it should be possible 
to implement uniform scaling around a given euclidean point in this way. 

With these remarks we close our discussion of the projective properties of 
bivectors and the induced null polarities. We'll meet the null system again in 
Section 6 since it is fundamental to understanding rigid body mechanics. 

Metric properties of bivectors. Write the bivector H as the sum of 
two simple bivectors H = Hoc + So: 

3oo := Poieoi + Po2eo2 + Poaeoa 
So := Pi2ei2 + Paieai + P23e23 

This is the unique decomposition of H as the sum of a line lying in the ideal 
plane (Hoc) and a euclidean part (Ho)- We sometimes write H = (Hoc; So). 
This decomposition is useful in characterizing the bivector. HqAHoo = <^=> 
Ho V Hoc = <^=> H is simple. We say a bivector is ideal if Ho = 0, otherwise 
it is euclidean. Ho is a line through the origin, whose direction is given by the 
ideal point N3 := epHo = P23E1 +P31E2 +P12E3. We call Ns the direction 
vector of the bivector. The following identities are left as exercises for the 
reader: 

Ns = eoH (23) 
Eo V eoHo = Ho (24) 

Ho is invariant under euclidean translations, Hqo is not. This deserves a 
closer look. First, consider the case of a line passing through the origin. Then 
H = Eg V N where N is an ideal point. Note that Hqo = in this case. Let 
T be the ideal point representing a translation vector. Then the image of H 
under this translation is: 

Hr = (Eo + T) V N (25) 
= H + TVN (26) 
= (TVN;Ho) (27) 

where we have used the fact that ideal points are invariant under translations, 
and that the join of two ideal points is an ideal line (T V N). This leads us 
to the following proposition: 

Theorem 2. A bivector of the form (T V epHo; Ho) where T is ideal, is a 
simple bivector passing through the point Eq + T, and every such bivector can 
be so represented. 
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Proof. Set N boHq in (27). Then Hoc = T V eoHo. The resulting bivector 
Hfx) + So is simple since V Hq = and passes through Eq + T: 

HV (Eo + T) = (TVeoHo + Ho)V (Eo + T) (28) 

= (TVeoHo) VEo + HoVT (29) 

= TV(EoVeoHo)-TVHo (30) 

= T V (Eo V eoHo - Ho) (31) 

= (32) 

Here we have used the fact that (T V eoHo) V T = 0, Ho V Eq = 0, anti- 
symmetry and distributivity of V, and, finally, (24). 



When H is not simple, we have a related result: 

Theorem 3. A bivector of the form (Hqo + T V epHo; Ho) where T is ideal, 
is the image of the bivector H = Hoc + Hq under the translation by T, and 
every such bivector can be so represented. 

Proof. We know from Thm. 2 that the result is true when Hqo = 0. But since 
Hoo is an ideal line, it is fixed by a euclidean translation. The result follows 
immediately. □ 

Thm. 2 can be so intepreted: the set of all lines sharing the same finite part 
Ho is the line bundle centered at the ideal point N^. This bundle is mapped 
to itself by a translation T, in such a way that two elements in the bundle 
differ by the ideal line bundle element T V egHo when one is the translated 
image of the other under T. This implies that a change of coordinate which 
moves Eq to Eg+T will produce such a result on the coordinates of bivectors. 
We'll return to this point below in the discussion of the axis of a bivector 
(Sect. 5.3), where we show how to choose a canonical representative from this 
line bundle to represent a general bivector H. . 



5.2 Enumeration of various products 

All the products described in 4.1 have counter-parts here, obtained by leaving 
points alone and replacing lines by planes. We leave it as an exercise to the 
reader to enumerate them. Here we focus on the task of enumerating the 
products that involve bivectors. For that purpose, we extend the definition 
of a,b, P, and Q to have an extra coordinate, and introduce two arbitrary 
bivectors, which may or may not be simple: H := poieoi + ... and $ := 

501601 + 

1. Inner product. H-$ = —{P12912+P31913+P23923) = — cos (a) where a is 
the angle between the direction vectors of the two bivectors (see Sect. 5.1 
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above). H • $ is a symmetric bilinear form on bivectors, called the Killing 
form. We sometimes write H • $ = (H, Note that just as in the 2D 
case, the ideal elements play no role in this inner product. This angle 
formula is only valid for euclidean bivectors. 

2. Norm. There are two cases: 

a. Euclidean bivectors. For euclidean H, define the norm |1H|| — H • H. 
Then norm 1; we call it a normalized euclidean bi vector. 

b. Ideal bivectors. As in Sect. 4.1, we get the desired norm on an ideal 
line H by joining the line with any euclidean point P, and taking the 
norm of the plane: ||H||oo — ||HVP||. We normalize ideal bivectors with 
respect to this norm. 

3. Distance. Verify that the euclidean distance of two normalized points P 
and Q is still given by ||P V Q|j, and the norm of an ideal point V (i. e., 
vector length) is given by ||V V P|j where P is any normalized euclidean 
point. 

4. Inverses. For euclidean H, define = — — —. Inverses are unique. 

5. HA# = HV$ = (H, $)pl is the Pliicker inner product times I. When 
both bivectors are simple, this is proportional to the euclidean distance 
between the two lines they represent (Exercise). 

6. Commutator. H x $ is a bivector which is in involution to both H and 
$ (Exercise). We'll meet this later in the discussion of mechanics (Sect. 6) 
as the Lie bracket. 

7. Null point, a A H, for simple H, is the intersection point of H with the 
plane a; in general it's the null point of the plane with respect H. 

8. Null plane. P V H, for simple H, is the joining plane of P and H; in 
general it's the null plane of the point with respect to H. 

9. P ■ H, for simple H, is a plane passing through P perpendicular to H. 

10. P X H, for simple H, is the normal direction to the plane through P and 

11. a • H, for simple H, is a plane containing H whose intersection with a is 
perpendicular to H. 

12. HI — (p23eoi +P3ieo2 +^12603)1 is the polar bivector of the bivector H. 
It is an ideal line which is orthogonal (in the elliptic metric of the ideal 
plane, see Sect. 3.1) to the direction vector of H. (Hoc; 3o)I — (SqI; 0). 



5.3 Dual Numbers 

A number of the form a + bl for a, & € R we call a dual number, after Study 
([Stu03]). Dual numbers are similar to complex numbers, except = rather 
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than P — —1. We'll need some results on dual numbers to calculate rotor 
logarithms below. 

• Dual numbers commute with other elements of the Clifford algebra. 

• Given a dual number z = a + 51, we say z is euclidean ii a ^ 0, otherwise 
z is ideal. 

• Conjugate. Define the conjugate z = a — bl. zz = a^. 

• Norm. Define the norm \\z\\ := \/z^ = a. 

• Inverse. For euclidean z, define the inverse z — The inverse is the 
unique dual number w such that zw = 1. 

• Square root. Given a euclidean dual number a + 61, define c = and 

d = — —. Then w := c + dl satisfies = z and we write w ^/z. 

Dual analysis. Just as one can extend real power series to complex power 
series with reliable convergence properties, power series with a dual variable 
have well-behaved convergence properties. See [Stu03] for a proof. In partic- 
ular, the power series for cos {x + yl) and sin {x + yl) have the same radii of 
convergence as their real counterparts. One can use the addition formulae for 
cos and sin to show that ([PWOl], pl54): 

cos {x + yT) = cos x — (yl) sinx 
sin {x + yV) = sin x + (yl) cos x 

The axis of a bivector. Working with euclidean bivectors is simplified 
by identifying a special line, the axis, the unique euclidean line in the linear 
span of Hoc and Ho- The axis is defined by = {a + 61) H for a dual 
number a + bl. In fact, one can easily check that the choice 

a:6=-2(H,HI)p:(H,H)p 

yields the desired simple bivector. We usually normalize so that — — 1. 
The axis appears later in the discussion of euclidean isometrics in Sect. 5.5, 
since most isometries are characterized by a unique invariant axis. 

5.3.1 Exercises 

1. Translate the products a • b, P • Q, a A P, a • P, P V Q, P x Q, al, PI from 
Sect. 4.1. [Hint: a typical 2D item can be translated to 3D by replacing lines 
by planes, and 2- vectors with 3- vectors. The grade of a product depends, as 
before, on the grades of the arguments. For example, a-P is the line passing 
through P perpendicular to the plane a.] Also, translate the definition of 
norm to points and planes in 3D. 

2. Relationships of lines. In this exercise, both H and $ are normed eu- 
clidean simple bivectors. 
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a. Distance between lines. For normalized H and HA^" = sin {a)ds^, 
where ds^ is the euclidean distance between the two lines, and a is the 
angle between their two direction vectors. [Hint: consider the tetrahe- 
dron spanned by unit vectors on H and €>.] 

b. Dual norm and dual angle. Define the dual norm ||H||d :— (H)o + 
(3)4. Then j|H$||rf = H and # intersect at right angles. In 
general, 

||H#||(i ~ ±(cos (a) — sin {a)ds^l 

where a is the angle between the direction vectors of H and This is 
called the dual angle of H and $ and measures both the angle between 
the directions and the distance between the lines. Give an algorithm to 
decide which choice of ± is correct. 

3. The axis. Let the axis of a non-simple H = Hoc + So be a^- Show that 

is the unique euclidean line such that Hqo is both the polar line of 
with respect to the euclidean metric (i.e.. Hoc = SI) and the conjugate 
line of Hix with respect to the null polarity on S. 

4. Orthogonal projection. It's interesting to investigate orthogonal pro- 
jectsions involving lines. Not only can one project points onto lines, but 
lines can also be projected onto points. To be more precise: One can project 
points onto spears, and axes onto bundles. Keep in mind in the following 
exercises that X"-'^ can in general be replaced by X without effecting the 
validity of the result, qua subspace. 

a. Projecting a point onto a line, and vice-versa. Show that for 
euclidean P and euclidean simple S, (P-S)S^^ is the point of S closest 
to P. What is ((P • S)S) V P? Show that (P • S)P-i is a fine parallel 
to S passing through P. 

b. Projecting a line onto a plane, and vice-versa.. Show that for 
euclidean a and simple S, (S • a)a^^ is the orthogonal projection of S 
onto a. Show that (a • S)a^^ represents a plane containing S, parallel 
to a. 

5. The common normal of two lines. Let S and $ be two simple eu- 
clidean bivectors. 

a. Show that := S x $ is a bivector which is in involution to both S 
and 

b. Show that the axis &x of is a line perpendicular to both S and 

c. Calculate the points where 02, intersects S and [Hint: consider where 
the plane spanned by P := BQ&r^ (the direction vector of &r^) and $ 
cuts S.] 

6. *Many of the products in Sect. 5.2 are described only for simple S. Can 
you provide an interpretation for non-simple S? 
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7. Assume H A 4' = for simple bivectors H and ^, with H 7^ Find the 
unique common point and unique common plane of H and [Hint: for a 
point P ^ H, ^ consider H A (* V P).] 

8. The pitch of a bivector. Let H be a euclidean bivector and define the 

pitch p to be the ratio: p :— —j — ' — ^r^. To account for ideal H, one defines 
p — 00. Show that 

a. For simple euclidean H, p = 0. 

b. p is a euclidean invariant. 

9. For two general bivectors H and use the ideal- finite decomposition 
H ~ Hoc + Ho and $ = $00 + *o (Sect. 5.1). Establish the following 
formulas: 

a. H A * = Hoo A *o + 3o A *oo- 

b. H X * = (Hoc X $0 + Ho X *oo; So X *o)- 

c. H • * = Ho • *o. 

d. H is simple <^=4> Hqo A Ho = 0. 

10. Prove the following identities involving involving general bivectors, and an 
arbitrary 3- vector P: 

(H X A H = (33) 

S((HI)A*) = S(HA(*I)) = H-* (34) 

H V (HAP) = S(HAH)P (35) 

HxP = (HVP)I (36) 

P X (H X P) =PHP + H (37) 

11. Incidence of point and line. Show that a euclidean point P lies on 
the simple bivector H <^=4> P V H = 0. Use (36) above to show this is 
equivalent to P x H = 0. Show that in both cases the formulas are valid 
also for non-simple H. (That is, for non-simple H, neither expression can 
vanish.) 



5.4 Reflections, Translations, Rotations, and ... 

The results of Sect. 4.2 can be carried over without significant change to 3D: 

1. For a 1-vector a, the sandwich operation P aPa is a euclidean refiection 
in the plane represented by a. 

2. For a pair of 1-vectors a and b such that g :— ab, P gPg is a euclidean 
isometry. There are two cases: 
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a. When (g)2 is euclidean, it's a rotation around the hne represented by 
(g)2 by twice the angle between the two planes. 

b. When (g)2 is an ideal line poieoi +Po2eo2 +Po3eo3, it's a translation by 
the vector {x,y,z) = 2(poi,Po2,Po3)- 

A rotor responsible for a translation (rotation) is called, as before, a translator 
{rotator). There are however other direct isometries in euclidean space besides 
these two types. 

Definition 5. A screw motion is a isomctry that can be factored as a rotation 
around a line H followed by a translation in the direction of H. H is called 
the axis of the screw motion. 

Like the linear line complex, a screw motion has no counterpart in 2D. In 

fact, (g)2 is a non-simple bivector <;=^ g is the rotor of a screw motion. To 
show this we need to extend 2D results on rotors. 



5.5 Rotors, Exponentials and Logarithms 

As in Sect. 4.3. the spin group Spin(3,0, 1) is defined to consist of all ele- 
ments g of the even subalgebra PCMs'o i) such that gg = 1. A group element 
is called a rotor. In this section we seek the logarithm of a rotor g. Things are 
complicated by the fact that the even subalgebra includes the pseudo-scalar 
I. Dual numbers help overcome this difiiculty. 
Write g = Sr + s^I + H. Then 

gg = + 2SrSdI - H2 = 1 

Suppose is real. Then H is simple, sj. = 0, and = — 1. If < 0, 
then find real A such that Hjv := AH satisfies = — 1, and evaluate the 
formal exponential e*^~ as before (Sect. 4.3) to yield: 

e*^'' = cos (i ) + sin (t ) H jv (38) 

We can use this formula to derive exponential and logarithmic forms for 
rotations as in the 2D case (Exercise). If = 0, H is ideal, and the rotor is 
a translator, similar to the 2D case (Exercise). This leaves the case Sd ^ 0. 
Let # = (a -I- 5I)H be the axis of H (see Sect. 5.3 above). Since the axis is 
euclidean, a ^0, and the inverse exists: c + dl := (a + bl)~^. 

(c + (n)* = H (39) 

Replace the real parameter t in the exponential with a dual parameter t + ul, 
replace H with substitute = —1, and apply the results above on dual 
analysis: 
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= COS (t) - u sin (<)I + (sin (<) + u cos (41) 

We seek values of t and u such that g equals the RHS of (41): 

g ^ Sr + Srfl + (c + dl)* = COS (t) - usm{t)I + {am{t) + ucos(i)I)* (42) 

We solve for t and u, doing our best to avoid numerical problems that might 
arise from cos (t) or sin (t) alone: 

t — tan~"'^(c, Sr) 

if |cos(f)| > |sin(t)| 



cos [t) 

— ^ otherwise 
. sm (t) 

Definition 6. Given a rotor g e 5*^(3, 0, 1) with non-simple bivector part, 
the bivector [t + mI)€> defined above is the logarithm of g. 

Theorem 4. Let {t + mI)4> be the logarithm of the rotor g (= + Sdl + 3) 
with Sd ^ 0. Then u ^ 0,t ^ and g represents a screw motion along the 
axis $ consisting of rotation by angle 2t and translation by distance 2u. 

Proof. $ commutes with H and with g (Exercise). Hence g€>g = $ is fixed 
by the sandwich. Write the sandwich operation on an arbitrary blade x as 
the composition of a translation followed by a rotation: 

gxg = e(*+"i)*xe-(*+"i)* 
= e**(e"^*xe-"^*)e-** 

This makes clear the decomposition into a translation through distance 2u 
(Exercise), followed by a rotation around ^ through an angle 2t. One sees 
that the translation and rotation commute by reversing the order. 

Note that calculating the logarithm of a rotor involves two separate nor- 
malizations. First, the rotor must be normalized to satisfy gg = 1. Then 
the axis of the bivector of g must also be extracted, a second normalization 
step. This axis, multiplied by an appropriate "dual angle" t + ul can then be 
exponentiated to reproduce the original rotor. 

Translations. Translations represent a degenerate case. The logarithm of 
a translator is not unique, since 



for any ideal bivector H. This is related to the fact that a translator has no 
well-defined axis, since the pencil used to define the axis (Sect. 5.3) degener- 
ates to the translator itself. But this degeneracy does not cause difficulties. 
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since the calculation of exponential and logarithms in this case is simpler 
than the general case. We choose the logarithm that goes through the origin 
of the coordinate system. 

We have succeeded in showing that the bivector of a rotor is non-simple 
if and only if the associated isometry is a nondegenerate screw motion. This 
result closes our discussion of P(IR3 q i). For a fuller discussion, see [Gunllb]. 

5.5.1 Exercises 

1. Handle the case Sd = (g = + H) from Sect. 5.5 to obtain exponential 
forms for rotations and translations. Define the corresponding logarithms. 
What does the case 3^ = 3^ = represent? 

2. Find the rotator corresponding to a rotation of — radians around the line 

through the origin and the point (1,1,1). [Answer: .5(1 + ei2 +e3i +e23).] 

3. From the proof of Thm. 4: 

a. Show that $ commutes with H and with g. 

b. Confirm the claim of the theorem that the translation moves points a 
distance 2u. [Hint: = —1 implies ||I$||oo = 1-] 



6 Case Study: rigid body motion 

The remainder of the article shows how to model euclidean rigid body motion 
using the Clifford algebra structures described above. It begins by showing 
how to use the Clifford algebra P(M3 q i) to represent euclidean motions and 
their derivatives. Dynamics is introduced with newtonian particles, which 
are collected to construct rigid bodies. The inertia tensor of a rigid body 
is derived as a positive definite quadratic form on the space of bivectors. 
Equations of motion in the force-free case are derived. In the following, we 
represent velocity states by O, momentum states by 11, and forces by A. 



6.1 Kinematics 

Definition 7. A euclidean motion is a path g : [0, 1] Spin(3, 0, 1) 
with g{0) = 1. 

Theorem 5. For a euclidean motion g, gg is a bivector. 

Proof, gg is in the even subalgebra. For a bivector X, X = —X; for scalars 
and pseudoscalars, X = X. Hence it suffices to show gg = — gg. 
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Fig. 6 For n = 2 the 

euclidean velocity state 
is a point, and it acts on 
points. The vector field 
2V X P in the neighbor- 
hood of V. 




gg = 1 



(gg) = 



gg + gg = 
gg + gg = 



gg = -gg 



□ 



Define := g(0); by the theorem, f2 is a bivector. We call ft a euclidean 
velocity state. For a point P, the motion g induces a path P(t), the orbit of 
the point P, given by P(t) = g{t)Pg{t). Taking derivatives of both sides and 
evaluating at i = yields: 



The last step follows from the definition of the commutator product of bivec- 
tors. In this formula we can think of P as a normalized euclidean point which 
is being acted upon by the euclidean motion g. From Sect. 5.2 we know that 
f2 X P is a ideal point, that is, a free vector. We sometimes use the alterna- 
tive form r2 X P = (r2 V P)I (Exercise). The vector field vanishes wherever 
Jl V P = 0. This only occurs if is a line and P lies on it. The picture is 
consistent with the knowledge, gained above, that in this case e*^ generates 
a rotation (or translation) with axis fl. Otherwise the motion is an instanta- 
neous screw motion around the axis of fl and no points remain fixed. Fig. 6 
shows how the vector field looks in case n = 2. It's easy to see that e*^ in 
this case yields a rotation around the point V. 



m 
p{t) 

P(0) 



s{t)Pg{t) + g(<)Pg(i) 
g(t)Pg(t) + g(i)Pi(i) 

r2P po 

2(n X p) 
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Null plane interpretation. In the formulation P = 2{fi V P)I, we rec- 
ognize the result as the polar point (with respect to the euclidean metric) 
of the null plane of P (with respect to Jl). See Fig. 6.1. Thus, the vector 
field can be considered as the composition of two simple polarities: first the 
null polarity on $7, then the metric polarity on the euclidean quadric. This 
leads to the somewhat surprising result that regardless of the metric used, 
the underlying null polarity remains the same. One could say, for a given 
point, its null plane provide a projective ground for kinematics, shared by all 
metrics; the individual metrics determine a different perpendicular direction 
to the plane, giving the direction which the point moves. This decomposition 
only makes itself felt in the 3D case. In 2D, the null polarity is degenerate: 
VVP is the joining line of V and P (a similar degeneracy occurs in 3D when 
n is simple). 



6.2 Dynamics 

With the introduction of forces, our discussion moves from the kinematic level 
to the dynamic one. We begin with a treatment of statics. We then introduce 
newtonian particles, build rigid bodies out of collection of such particles, and 
state and solve the equations of motions for these rigid bodies. 

See Appendix 3 for a detailed account of how 2D statics are handled in 
this model. 

3D Statics. Traditional statics represents 3D a single force as a pair 
of 3- vectors {V,M), where V = {vx,Vy,Vz) is the direction vector, and M — 
{mx,rny,mz) is the moment with respect to the origin (see [Fea07], Ch. 2). 
The resultant of a system of forces Fi is defined to be 



^F, = (^y„^A/,)=: {V,M) 

i i i 
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The forces are in equiUbrium <;=4> V — M = 0.V" = and M ^ 

the resultant force is a force couple. Otherwise the vectors V and M are 

orthogonal ■<=^ the system represents a single force. 

If P is a normalized point on the line carrying the force, define H{F) :~ 
P V i{V). We call H{F) the homogeneous form of the force, and verify that: 

H{F) = TO^Eoi + m.yEo2 + m^Eos + WzEi2 + WyEgi + w^Eas 

If F is the resultant of a force system {F,}, then H{F) = H{Fi). Hence, 
a system of forces {Fi] is the null force H{Fi) — 0. Furthermore, 

H{F) is an ideal line <^=^ the system of forces reduces to a force-couple, 
and H{F) is a simple euclidean bivector <;=^ F represents a single force. 
Notice that the intensity of a bivector is significant, since it is proportional 
to the strength of the corresponding force. For this reason we sometimes say 
forces are represented by weighted bivectors. 



6.2.1 Newtonian particles 

The basic object of Newtonian mechanics is a particle P with mass m located 
at the point represented by a trivector R. Stated in the language of ordinary 
euclidean vectors, Newton's law asserts that the force F acting on P is: 
F = mR. 

Definition 8. The spear of the particle is A := R V R. 
Definition 9. The momentum state of the particle is II :— mA. 
Definition 10. The velocity state of the particle is F := AI. 
Definition 11. The kinetic energy E of the particle is 

TTi • m 1 

E := -IIRIIL = -yA . A = --S(r A U) (43) 

Remarks. Since we can assume R is normalized, R is an ideal point. 11 is 
a weighted bivector whose weight is proportional to the mass and the velocity 
of the particle. T is ideal, corresponding to the fact that the particle's motion 
is translatory. Up to the factor m, T is the polar line of II with respect 
to the euclidean metric. It is straightforward to verify that the linear and 
angular momentum of the particle appear as Ho and IIoo, resp., and that 
the definition of kinetic energy agrees with the traditional one (Exercise). 
The second and third equalities in (43) are also left as an exercise. 

We consider only force-free systems. 

Theorem 6. If F — then A, II, T, and E are conserved quantities. 



Proof. F — implies R = 0. Then: 
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• A = (RVR + RVR)=0 

• n = mA = 
. f = (A)I ^ 

• £; = ^(S(fAn) + s(rAri)) = 0. 

Inertia tensor of a particle. Assume the particle is "governed by" a 



euclidean motion g witlr associated euclidean velocity state fi := g(0). Then 
n, r, and E depend on fl as follows: 

R = 2(J7 X R) (44) 

n = 2m(RV (rj X R)) (45) 

r = 2(RV (J7 X R))I (46) 

= 2(R x{nx R)) (47) 

£; = -2mS(rAn) (48) 

= -2toS((R X (J7 X R)) A (RV X R))) (49) 

= -S{nAU) (50) 



The step from (46) to (47) follows from (37). The step from (48) to (50) is 
equivalent to the assertion that 2(11 A 11) = F A 11. From 11 ~ mR V R it 
is enough to show 2(r2 V R) = F V R. Since both sides of the equation are 
planes passing through R, it only remains to show that the planes have the 
same normal vectors. This is equivalent to 

Theorem 7. 2(f2 x R) = F x R. 

Proof. 

FxR = 2(Rx (nxR) xR) (51) 

= 2(R(fl X R)R) (52) 

= (R(rjR - Rf2)R) (53) 

= (R0R2 - R^OR) (54) 

= (-Rr2 + nR) (55) 

= 2{n X R) (56) 

Here we have used the fact that R • (Jl x R) = 0, the definition of x , the fact 
that for euclidean points R^ = — 1, and finally the definition of x a second 
time. □ 

Define Then the theorem yields immediately two corollaries: 

Corollary 1. The difference H := 2fl — T is a simple bivector incident with 
R. 



Proof. The theorem implies H x R = 0. Since this is the normal direction of 
the plane H V R, and H is euclidean, this implies H V R = 0. The null plane 
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of a point, however, can only vanish if the bivector is simple and is incident 
with the point. □ 

Corollary 2. H is the conjugate line ofT with respect to the null polarity fl. 

Proof. This follows from the observation that the conjugate of a line $ with 
respect to a non-simple bivector f2 is a line lying in the pencil spanned by # 
and fl. This condition is clearly satisfied by both H and T. Since there are 
at most two such lines in the pencil, the proof is complete. □ 

We return to the theme of Newtonian particles below in Sect. 6.2.2. 
Define a real-valued bilinear operator A on pairs of bivectors: 

777 

A{n,a) — S(((RV2(OxR))I)A(RV2(HxR))) (57) 
ni 

= — (RV2(f2 X R)) • (RV2(H X R)) (58) 

where the step from (57) to (58) can be deduced from Sect. 5.2. (58) shows 
that A is symmetric since • on bivectors is symmetric: A • A = A • A. We call 
A the inertia tensor of the particle, since E ~ A(r2, O) = —8(11 A 11). We'll 
construct the inertia tensor of a rigid body out of the inertia tensors of its 
particles below. We overload the operator and write 11 = A(r2) to indicate 
the polar relationship between 11 and CI. 

6.2.2 Rigid body motion 

Begin with a finite set of mass points Pf, for each derive the velocity state 
Ti, the momentum state 11^, and the inertia tensor A.^.^^ Such a collection of 
mass points is called a rigid body when the euclidean distance between each 
pair of points is constant. 

Extend the momenta and energy to the collection of particles by summa- 
tion: 

n:=^n, = ^A,(rj) (59a) 
£; := ^B,; = ^ Ai(J7,r2) (59b) 

Since for each single particle these quantities are conserved when F = 0, 
this is also the case for the aggregate 11 and E defined here. 
We introduce the inertia tensor A for the body: 

Definition 12. A := X] A^. 



We restrict ourselves to the case of a finite set of mass points, since extending 
this treatment to a continuous mass distribution presents no significant technical 
problems; summations have to be replaced by integrals. 
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Then 11 = A(r2) and E = A{fl,fl), neither formula requires a summation 
over the particles: the shape of the rigid body has been encoded into A. We 
sometimes use the identity 

A(f2i, fla) = -S(rJi A A(ri2)) (Exercise) (60) 

which is a consequence that the individual inertia tensors for each particle 
exhibit this property. 

One can proceed traditionally and diagonalize the inertia tensor by finding 
the center of mass and moments of inertia (see [Arn78]). Due to space con- 
straints we omit the details. Instead, we sketch how to integrate the inertia 
tensor more tightly into the Clifford algebra framework. 

Clifford algebra for inertia tensor. We define a Clifford algebra Ca 
based on P(/\^]R^*) by attaching the positive definite quadratic form A as 
the inner product. We denote the pseudoscalar of this alternative Clifford 
algebra by Ia, and inner product of bivectors by (,)a- We use the same 
symbols to denote bivectors in W* as I-vectors in Ca. Bivectors in W are 
represented by 5-vectors in Ca- Multiplication by I a swaps I-vectors and 5- 
vectors in Ca; we use J (lifted to Ca) to convert 5-vectors back to I-vectors 
as needed. The following theorem, which we present without proof, shows 
how to obtain 11 directly from I a in this context: 

Theorem 8. Given a rigid body with inertia tensor A and velocity state rj, 
the momentum state 11 = A(J1) = J(rilA). 

Conversely, given a momentum state II, we can manipulate the formula in 
the theorem to deduce: 

o = A-i(n) = (j(n)ix^) 

In the sequel we denote the polarity on the inertia tensor by A(ll) and 
A~^(n), leaving open whether the Clifford algebra approach indicated here 
is followed. 

Newtonian particles, revisited. Now that we have derived the inertia 
tensor for a euclidean rigid body, it is instructive to return to consider the 
formulation of euclidean particles above (Sect. 6.2.1). We can see that in this 
formulation, particles exhibit properties usually associated to rigid bodies. 

• E = -^S{T A n): The kinetic energy is the result of a dual pairing be- 
tween the particle's velocity state and its momentum state, considered as 
bivectors. 

• r = ^ni: The dual pairing is given by the polarity on the euclidean 
metric quadric, scaled by ^. This pairing is degenerate and only goes in 
one direction: from the momentum state to produce the velocity state. 



It remains to be seen if this approach represents an improvement over the Unear 
algebra approach which could also be maintained in this setting. 
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• E — S{fl A n): the same energy is obtained by using twice the global 
velocity state in place of the particle's velocity state. This follows from 
Thm. 7. 



6.2.3 Exercises 

1. Verify that the linear and angular momentum of a particle appear as IIo 
and IIoo, resp. 

2. Verify the equalities in (43). 



6.2.4 The Euler equations for rigid body motion 

In the absence of external forces, the motion of a rigid body is completely 
determined by its momentary velocity state or momentum state at a given 
moment. How can one compute this motion? First we need a few facts about 
coordinate systems. 

Coordinate systems. Up til now, we have been considering the behavior 
of the system at a single, arbitrary moment of time. But if we want to follow a 
motion over time, then there will be two natural coordinate systems. One, the 
body coordinate system, is fixed to the body and moves with it as the body 
moves through space. The other, usually called the space coordinate system, 
is the coordinate system of an unmoving observer. Once the motion starts, 
these two coordinate systems diverge. The following discussion assumes we 
observe a system as it evolves in time. All quantities are then potentially 
time dependent; instead of writing g{t), we continue to write g and trust the 
reader to bear in mind the time-dependence. 

We use the subscripts Xg and Xc^^ to distinguish whether the quantity X 
belongs to the space or the body coordinate system. The conservation laws of 
the previous section are generally valid only in the space coordinate system, 
for example. Tig = 0. On the other hand, the inertia tensor will be constant 
only with respect to the body coordinate system, so, lie = A(rJc). When we 
consider a euclidean motion g as being applied to the body, then the relation 
between body and space coordinate systems for any element X S P(M3 q i), 
with respect to a motion g, is given by the sandwich operator: 

X, = gX,g 

Definition 13. The velocity in the body flc gg, and the velocity in space 
fls := g^cg. 

Definition 14. The momentum in the body Tl^ := A(Jlc), and the momen- 
tum in space Us := gllcg. 



From corpus, Latin for body. 
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We derive a general result for a time-dependent element (of arbitrary 
grade) in these two coordinate systems: 

Theorem 9. For time-varying X G P(]R3Q]^) subject to the motion g with 
velocity in the body Q,c, 

X, = g(X, + 2(f2, X X,))g 

Proof. 

Xs = gX,g + gX,g + gX,i 

= g(ggXc +Xc+ Xcgg)g 

= g(r2cXc + Xc + Xcf2c)g 

= g(X, + J^eX, - X,n,)g 

= g(X, + 2(0, xX,))g 

The next-to-last equality follows from the fact that for bivectors, Q = — f2; 
the last equality is the definition of the commutator product. □ 

We'll be interested in the case Xc is a bivector. In this case, Xc and ric 
can be considered as Lie algebra elements, and 2{flc x Xc) is called the Lie 
bracket. It expresses the change in one (X) due to an instantaneous euclidean 
motion represented by the other ($7). 

6.2.5 Solving for the motion 

Since fic = gg, g = S^c, a first-order ODE. If we had a way of solving for 
r2c: we could solve for g. If we had a way of solving for 11^, we could apply 
Theorem 8 to solve for ftc- So, how to solve for He? 

We apply the corollary to the case of force-free motion. Then 11^ = 0: the 
momentum of the rigid body in space is constant. By Theorem 9, 

= tls=g{Uc + 2{n,xU,))g (61) 

The only way the RHS can be identically zero is that the expression within 
the parentheses is also identically zero, implying: 

lie = X r2c 

Use the inertia tensor to convert velocity to momentum yields a differential 
equation purely in terms of the momentum: 

lie = X A-i(nc) 

One can also express this ODE in terms of the velocity state alone: 
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lie = A-\u,) = 2A-i(n, X A-i(nc)) 

= 2A"i(A(Oc) X flc) 



When the inner product is written out in components, one arrives at the 
well-known Euler equations for the motion ([Arn78], p. 143). 

The complete set of equations for the motion g are given by the pair of 
first order ODE's: 



where flc — A^^(nc). When written out in full, this gives a set of 14 first- 
order linear ODE's. The solution space is 12 dimensions; the extra dimensions 
corresponds to the normalization gg = 1. At this point the solution continues 
as in the traditional approach, using standard ODE solvers. Our experience 
is that the cost of evaluating the Equations (62) is no more expensive than 
traditional methods. 



6.2.6 External forces 

The external force F acting on the rigid body is the sum of the external forces 
Fi acting on the individual mass points. In traditional notation, Fi = rriir. 
What is the homogeneous form for Fi? In analogy to the velocity state of 
the particle, we define the acceleration spear of the particle :— V Ri, 
and the force state of the particle := rriiTi. Then the total force A is the 
sum: 



The following theorem will be needed below in the discussion of work 
(Sect. 6.2.7): 

Theorem 10. A = H 

Proof. Take the derivative of (59a): 



g = gf^c 

ric = 2nc X vtc 



(62) 
(63) 




ETOi(R, VR,) 
A 



□ 
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Theorem 10 applied to (61) yields Euler equations for motion with external 
forces:: 

A, = = g(ri, + 2{n, x njji 
gA,g = lie + 2{n, X 
A, = u, + 2{n, X nj 
ric - A, + 2(nc X n^) 

Note that the forces have to converted from world to body coordinate systems. 
6.2.7 Work 

As a final example of the projective approach, we discuss the concept of work. 
Recall that 11 — A(r2), so 11 = A(r2), and the definition of kinetic energy 
for a rigid body : E = ^A(0, n). 

Theorem 11. E = -S{n A A) 

Proof. 

Til 

= mA{n,n) 

= -S(J7 A ?7iA(0)) = -S(0 A mil) 
= -S(f2A A) 

where we apply Leibniz rule, the symmetry of A, (60) and finally Thm. 10. 
□ 

In words: the rate of change of the kinetic energy is equal to the signed 
magnitude of the outer product of force and velocity. This is noteworthy in 
that it does not involve the metric directly. 

E is sometimes called the power. The work done by the force between time 
to and t is the defined to be the integral: 

w{t) = E{t) - E{to) = / Eds 

Jta 

= [ -S(f2AA)ds 

J to 

The integrand depends on the incidence properties of the force and the veloc- 
ity, as a line/point pair. If the two elements are incident, then f A = and 
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there is no work done; the further away ft lies from A, (m *B!) the greater 
the power and hence the work done. One can be more precise: 

Theorem 12. Suppose A. is a single force, Q, is a rotator, d is the euclidean 
distance between the lines A and il, and a is the angle between the two 
direction vectors. Then S(0 A A) — — (isinQ;||ll||||A||. 

Proof. Left as an exercise. □ 

In words: the rate of change of the kinetic energy is proportional to the 
intensity of the force, the intensity of the velocity, and the euclidean distance 
between the force line and the velocity line. 

Example. Imagine an ice skater who moving along the surface of a frozen 
lake with negligible friction; the single force is given by gravity. Assuming 
gravity is in the negative z-direction acting on the skater located at the 
origin, then A = gei2 (corresponding to the intersection line of the planes 
X = and y — with weight gravitation constant g). Consider two possible 
motions of the skater: 

• The motion of the skater is a translation in the x-direction given by an 
ideal line of the form 17 = deoi, d < 0. AAJ7 = 0, so no work is required 
for the skater to skate! 

• The skater spins around a fixed point. Then the velocity state relative 
to the natural diagonalized form of the inertia tensor has null ideal part 
rJoo — and the corresponding momentum state 11 = A(f2) has null 
euclidean part TIq = 0: it's a momentum couple: a momentum carried by 
an ideal line! 

• As the skater spins, she stretches her arms out, then pulls her arms close to 
her body. This latter movement decreases the entries in the inertia tensor 
A, increasing the entries in A~^; since fl = A~^(n), her velocity increases 
proportionally in intensity: she spins faster. 

One can see from this example the advantages of the projective approach: it 
unifies calculations, and handles exceptional cases, such as the translations 
and couples in the above example, at no extra cost. 

Comparison. The projective Clifford algebra approach outlined here ex- 
hibits several advantages over other approaches to rigid body motion. The 
representation of kinematic and dynamic states as bivectors rather than as 
pairs of ordinary 3- vectors (linear and angular velocity/momentum/etc.) pro- 
vides a framework free of the special cases which characterize the split ap- 
proach (for example, translations are rotations around an ideal line, a force 
couple is a single force carried by an ideal line, etc.). The Clifford algebra 
product avoids cumbersome matrix formulations and, as seen in Thm. 9, 
yields formulations which are valid for points, lines, and planes uniformly. 
The inertia tensor of a rigid body can be represented as a separate (posi- 
tive definite) Clifford algebra defined on the space of bivectors. Finally, the 
treatment of Newtonian particles reveals an underlying velocity-momentum 
polarity in bivector space analogous to that of rigid bodies. 
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7 Extension to elliptic and hyperbolic metrics 



The theoretical basis for modeling rigid body motion in noneuclidean spaces 
using linear line complexes was worked out in the 19*'* century and the earlier 
part of the 20*^* century ([Lin73], [Hea84], [dF02], [Bla42]). The treatment 
given here for euclidean geometry and physics can be extended to elliptic and 
hyperbolic space by changing the defining signature to be (4, 0, 0) for elliptic 
and (3, 1, 0) for hyperbolic space. Surprisingly little else needs to be changed. 
Fig. 7 shows a simulation of rigid body motion in hyperbolic space using the 
Clifford algebra Cl{3, 1,0). A detailed description of the phenomenology lies 
outside the scope of this article. 

Fig. 8 Rigid body mo- 
tion in hyperbolic space. 
Note that the motion 
is bounded but has no 
fixed point, instead seems 
to fill out a solid torus. 
The momentum state is 
a hyper-ideal line which 
lies completely outside 
hyperbolic space. 



In this way one obtains a setting in which euclidean geometry is positioned 
as a limiting case of geometries with nondegenerate metric; the connections 
to its neighbors sheds new light on euclidean geometry and physics. Since 
Riemann, we know that the metric of space can only be learned through 
empirical research. The difference to nondegenerate metric spaces is partic- 
ularly interesting in the case of rigid body motion. Here one becomes aware 
of peculiarities of euclidean space. 

For example, one considers a rigid body as composed of points. But in a 
nondegenerate metric, every point has a unique polar plane which kinemati- 
cally and dynamically cannot be separated from it - moves along with it, as it 
were. The one-sided focus on points has no validity here. One of Lindemann's 
important contributions was, in this setting, to state the underlying dynam- 
ics of a rigid body, not in terms of the instantaneous lines of translation of 
points, but in terms of the instantaneous axes of rotation of planes. Another 
difference: the concept of center of mass is specific to euclidean space; in el- 
liptic space a rigid body has 4 centers. But space constraints prevent further 
discussion of this interesting theme. 
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8 Guide to the Literature 

The work presented here, in its conceptual basis, is derived from [Stu03]. This 
seminal book worked out in impressive detail the structure which in modern 
form appears as the even subalgebra P(]R3"q ^). Study avoided using the term 
quaternion; the structure he described has nonetheless become known as the 
dual quaternions. His parameter e maps to the pseudoscalar I. A full descrip- 
tion of the correspondence between the two systems lies outside the scope 
of this work, nor is the full scope of [Stu03] reflected in material presented 
here.[Zie85] gives an excellent survey of the historical development that led 
up to Study's work. Mobius, Pliicker, Hamilton and Klein were Study's most 
important predecessors; he and Ball ([BalOO]) had a relationship based on mu- 
tual appreciation. Weiss, a student of Study's, wrote a concise introduction 
[Wei35] to Study's investigations which can be recommended. For beginners, 
[Bla54] provides a simpler introduction to Study's approach. 

Study's contribution in the direction of mechanics was developed further 
in [vM24] and [Bla42]. The former was hampered by the awkward matrix 
notation required for stating transformation laws. These, on the other hand, 
are "built in" to the Clifford algebra approach and simplify the approach 
considerably. [Bla42] concentrated on the noneuclidean case. 

The modern legacy of this work is varied. Some modern literature, such as 
[Fea07] , use spatial vectors to model rigid body motion; these are 6D vectors 
equivalent to our bivectors, but developed within a linear algebra framework 
reminiscent of [vM24]. While the dual quaternions have developed a follow- 
ing in robotics and other applied areas ([McC90]), their embedding as the 
even subalgebra of P(M5.o.i) ^^.s not received widespread acceptance. Mod- 
ern sources using P(]R3g ^) include [SelOO] and [Sel05]. Much contemporary 
work which applies Clifford algebra methods to physics and engineering uses 
the conformal model ([DL03], [Per09], [HeslO]) rather than the homogeneous 
model presented here. Our handling of rigid body motion owes much to the 
spirit of [Arn78]. 

Study himself appears to be aware of the possibility of extending his work 
within a more comprehensive algebraic structure. He remarked at the end of 
[Stu03] (p. 595, translated by the author): 

The elementary geometric theory, that hovers thus before us, will surpass the 
construction possibilities of the quaternions, to the same degree that the ge- 
ometry of dynanien surpasses the addition of vectors. The accessory analytic 
machinery will consist of a system of compound quantities, with eight, or better 
yet, with sixteen units. {Study's italics!} 

It seems likely that P{M.l q i) is in fact the 16-dimensional algebraic realiza- 
tions of Study's prophetic inkling. 
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9.1 Loose ends 

There are some loose ends in the formulation presented here, that would be 
good to tidy up: 

• Orientation. The derivation of the reflection in Sect. 4.2 shows that, 
before dehomogenizing, the reflection actually flips the sign of the trans- 
formed coordinates. In particular the homogeneous coordinate is -1, not 1. 
One needs to first establish that this phenomenon is generally present in 
indirect versors (generally, odd-grade versors) and also for other types of 
operands (sandwich filling), not just points. If it is generally true, it would 
be interesting, since such a reversal of sign can be interpreted as a change 
of orientation due to the reflection. That would imply that dehomogenizing 
should be applied selectively, in order to retain this information. For back- 
ground on this theme see [Sto92]. Previous work on this theme ([Pap96]) 
used metric polarity to implement duality instead of the projective ap- 
proach used here, but presumably the results obtained do not depend on 
this fact. 

• Software implementations. The past decade has seen a explosion of 
software which implements geometric algebras of various sorts. Such soft- 
ware is generally of two kinds: implementations in so-called general pur- 
pose languages such as C-I--I- or Java, and those based on symbolic algebra 
systems such as Mathematica or Maple. Examples of the former include 
[Fon], [Per], and [Hil]; and of the latter, [JLAV]. These generally have not 
been written with degenerate metrics in mind. Consequently, either there 
is no support for degenerate metrics, or it hasn't been debugged. As the 
discussion has made clear, we recommend that the implementation of what 
is called duality to be based on projective principles, not on a metric po- 
larity (Appendix 1). This typically should not present a serious problem 
for software systems, since true duality can in fact be simulated using the 
metric polarity. One needs also to establish to what extent the assump- 
tion of a non-degenerate metric has been "hard-wired" into the code, for 
example, in the existence of I^^. The author adapted [JLAV] to support 
the degenerate metric used here; these changes were not difficult to imple- 
ment. In correspondence with the authors of the package, it appears these 
changes will be integrated in the next release of the package. 



9.2 Visualization 

One essential task is the visualization of non-simple bivectors. As the above 
hopefully makes clear, bivectors are the essential ingredient to almost every 
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non-trivial phenomena in 3D kinematics and dynamics. Visualization tools 
for non-simple bivectors are however not well-developed. Naively drawing all 
lines of the linear line complex associated to the bivector is not practical. If 
one takes a metric as given, then one can draw the axis of the bivector, along 
with some extra information to identify which element of the bivector pencil 
is intended. One also need to consider ways to distinguish between axes and 
spears visually. 



10 Evaluation and comparison 

We collect and present observations on the use of the homogeneous model 
of euclidean geometry and its Clifford algebraic implementation, which arose 
during the course of the work. 

• Fits the phenomena. The mathematical theory of rigid body motion 
maps precisely onto the projective model presented here. Indeed, a series 
of renowned mathematicians in the 19*'' century (Pliicker, Klein, Clifford, 
Study, and others), in their research on rigid body motion, laid the foun- 
dations of the mathematical theory which later developed into modern 
Clifford algebras ([Kle72], [Stu03]). 

• Duality. Our implementation provides duality via the projective founda- 
tion of the homogeneous model, rather than attempting to obtain it using 
a metric. Appendix 1 makes clear why this projectively valid method is 
often confused with a metric polarity. 

• Ideal elements play an important role in the homogeneous model of 
euclidean space. In incidence calculations involving parallel elements, they 
provide correct answers. Free vectors are represented as ideal points. An 
ideal line represents, in statics, a force couple; in kinematics, the axis of a 
euclidean translation, and plays an important role throughout the theory 
presented here. 

• Efficient. Before optimizing, a multiplication in the homogeneous model 
is 4 times cheaper than in the conformal model. Integration of equations of 
motion is easier in the homogeneous model since there are fewer directions 
to wander off the integral surface. The only possibility for "wandering off" 
when solving the equations of rigid body motion is normalizing the motion 
g so that gg = 1. 

• Intuitive. The geometric representation in the homogeneous model agrees 
with the naive sense perceptions of the human being. Straight is straight. 

• Insider's view of non-euclidean spaces is naturally computed. This im- 
portant visualization feature of the projective model lies outside the scope 
of this article; a detailed discussion can be found in [GunlO]. 

• Degenerate metric. Rather than considering this a disadvantage, our 
experience leads us to appreciate the benefits it brings. This is the metric 
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of euclidean space, and it provides correct formulas and results for eu- 
clidean geometry. While the non-invertiblity of I is unfamiliar, it presents 
no difficulties in the calculations and procedures represented here. 



11 Conclusion 

We have established that P(IR* q^) is a model for euclidean geometry. By 
using a mixture of projective, ideal and properly euclidean elements, we have 
avoided the problems traditionally associated with degenerate metrics. The 
result provides a complete and compact representation of euclidean geometry, 
kinematics and rigid body dynamics. We hope that this work will stimulate 
others working in these fields, to consider the homogeneous model as a prac- 
tical solution to their problem domain, and to deepen and extend these initial 
results. 



Appendix 1: J, metric polarity, and the regressive 
product 

Because of its importance in our approach to the homogeneous model, we 
provide here an more stringent, dimension-independent treatment of the map 
J and its close connection to the polarity 11 on the elliptic metric, concluding 
with reasons to prefer the use of J to that of 11 for implementing duality in 
the Clifford algebra setting. 

Canonical basis. A subset S = {ii,i2, ■■■ik} oi N = {1,2, ...,n} is called 
a canonical "k-tuple oi N if ii < 12 < ... < ik- For each canonical fc-tuple of TV, 
define to be the canonical (n — fc)-tuple consisting of the elements N\S. 
For each unique pair {5, S*-*-}, swap a pair of elements of 5'"'" if necessary 
so that the concatenation 5*5'-'-, as a permutation P of TV, is even. Call the 
collection of the resulting sets S. For each 5 G 6, define e'^ = e'i...e''=. We 
call the resulting set {e"^} the canonical basis for W generated by {e*}. 

Case 1: Equip W = P(/\(M")) with the euclidean metric to form the 
Clifford algebra P(IR„ 0.0) with pseudoscalar / = e^. Then, by construction, 
e^I = e*^"^^) (remember (e'^)~^ = e'^).This is the polarity 11 : W ^ W on 
the elliptic metric quadric (see Sect. 3.1). 

Case 2: Consider W* , the dual algebra to W. Choose a basis {ei, 62, ...e,i} 
for W* so that represents the same oriented subspace represented by the 
(n — l)-vector e*^* of W. Construct the canonical basis (as above) of W* 
generated by the basis {e^}. Then define a map J : W ^ W* by J(e^) ~ es± 
and extend by linearity. J is an "identity" map on the subspace structure of 
MP": it maps a fc-blade B € W to the (n — fc)-blade € W* which represents 
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the same geometric entity as B does in MP". Proof: By construction, 
represents the join of the 1-vectors e^^, {ij £ S) in W. This is however the 
same subspace as the meet of the n~ k basis 1-vectors e^^. , {ij e S-^) of W* , 
since contains e-' exactly when j i- 

Conclusion.: Both J and 77 represent vahd grade-reversing involutive 
isomorphisms. The only difference is the target space: 77 : P(VF) — P(W), 
while J : P{W) -> P{W*). 

The regressive product via a metric. Given the point-based exterior 
algebra with outer product A A B, (representing the join of subspaces A 
and B), the meet operator AV B (also known as the regressive product) is 
often defined n{n{A) A 77(B)) [HS87]. That is, the euchdean metric (any 
nondegenerate metric suffices) is introduced in order to provide a solution to 
a projective (incidence) problem. By the above, the same result is also given 
via J(J(A) A J (7?). (Here, the A denotes the outer product in W*). 

The Hodge ★ operator. An equivalent method for producing dual coor- 
dinates is described in [PWOl], p. 150. The * operator is presented as a way 
of generating dual coordinates, which is an apt description of the J operator 
also. One can also define the regressive product by -k{-kA A ★7?). Formally, 
however, since ★ is a map from to 14^, it is identical to the metric polarity 
77. 

Comparison. The two methods yield the same result, but they have very 
different conceptual foundations. As pointed out in Sect. 2.2, the meet and 
join operators live the exterior algebra W and its dual. J provides the bridge 
between these two projective algebras, hence provides a projective explana- 
tion for what is a projective operation. A related advantage of J is that it 
also is useful when used "alone" - whereas 77 (and also -k) is only valid when 
it appears in the second power, and hence disappears. For these reasons, we 
propose a differentiation of the terminology to reflect this mathematical dif- 
ferentiation. Instead of referring to multiplication by I (or I""'^) as the duality 
operator, we propose it should be called the metric polarity operator This 
is consistent with the mathematical literature. The purely projective term 
duality would be reserved for the J operator. 



Appendix 2: The euclidean distance function as a limit 

Rewrite (10) and manipulate: 

cos^ (d,)(X,X),(Y, Y), = (X, Y), (64) 
cos2 (d,)(e + x\){e + yf) = (e + {xm)? (65) 
cos^ (d,)(e' + e((yi)' + (xi)^) + i^mf) = + HxiVi) + {xnjif (66) 

Now consider the limit as e oo. It's clear from (10) that lime^oo cosd^ — 1. 
So we can replace cosde^ by (1 — c?^), and simplify the resulting expressions: 
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(1 - dl){e^ + e{yl + xl) + {xmf ) = + H^iVi) + {xmf (67) 
e{yl + xl- 2xm) = dl{e^ + e{yj + xj) + {x^yif) (68) 

(xi - y,f = dl{e + {yl + xl) + ^(xij/i)^) (69) 

Notice that the LHS is exactly the desired euchdean metric. In order to 
make the hmit converge to this, we define a new distance function :~ \J tdf_- 
This is exactly the scaling needed to prevent the distance from going to zero 
in the limit. One obtains: 

= ixi - yif (71) 



Appendix 3: 2D Statics 

We include here a detailed account of euclidean planar statics, to show how 
the classical separation of forces into a linear and an angular part is overcome 
through the use of the projective model. 

Statics is the study of systems of rigid bodies in equilibrium, or close to 
equilibrium.^^ This discussion is restricted to systems of forces acting on a 
single rigid body. In two dimensional statics, forces are characterized by the 
following geometrical properties: 

1. A single force can be represented as a vector free to move along a given 
lin. One speaks of the vector associated to the force, and of the carrying 
line of the force. 

2. Two forces acting on lines which intersect, are equivalent to a single force 
acting on a line passing through the intersection of the two lines, whose 
vector is equal to the vector sum of the vectors of the two original forces. 
See Figure 9. 

3. Two forces acting along the same line are equivalent to a single force acting 
on this line, whose vector is the vector sum of the two vectors. 

If the vector of a force vanishes, we say it is the null force. 

We recall some standard results of traditional treatments of 2D statics. 
Suppose V = {vx,Vy) is the vector of a force, P = {px,Py) is a point on the 
carrying line. 

Definition 15. The moment of a force F with respect to the origin is M :— 

Px^y Py^x- 



Dynamics on the other hand is concerned with an analysis of the motion which 
results when a system is not in equilibrium. 
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The main result of planar statics takes the following form: 

Theorem 13. A system of forces Fi is equivalent to the null force <f=^ 

i i 

where V; is the vector of Fi and Mi is the moment of Fi with respect to the 
origin. 

Let's find a homogeneous version of the above! Given a force F , with P 
and Q two normalized points on the line carrying the force, such that the 
force vector is given by Q — P = {v^, Vy). Define the 1-vector H{F) ;= P V Q. 

Theorem 14. A system of forces Fi is equivalent to the null force 
j:.HiF,)=0. 

Proof. For normalized points P and Q, P := Eq + PxEi + Pj,E2 and Q := 
Eo + feEi+(7.yE2, 

P V Q (j)a:qy - Pyqx)eo - [qy - Py)ei + (fe - Px)e2 (72) 

Note first that the coordinates of the force vector occur (rotated as it were 
by 90 degrees) as the coefficients of ei and e2. A simple computation shows 
that the coefficient of eo is the moment of the force m: 

PxQy - PyQx = Px{Py + Vy) - Py{px + Vx) 
= PxVy - PyVx 

Hence P V Q = meg + Vyei — Vxe2- The summand in this theorem vanishes 
«^=4> both summands in the previous one vanish. □ 

In brief: the blade determined by the force vector's head and tail, contains 
both force vector and moment with respect to the origin. This motivates the 
following definition: 

Definition 16. Let be a force with vector {vx,Vy) and moment to. Then 
H{F) — meo + VyBi — Vxe2 is called the homogeneous form of F. 
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Remarks. As the theorem shows us, the homogeneous form replaces two 
equations with a single one. Furthermore, it provides a unified form which 
subsumes a sticky exception in euclidean statics, the so-called force couple^ 
which arises as the resultant of two forces with equal-and-opposite force vec- 
tors carried on parallel lines. There is no single force equivalent to a force 
couple, so all results in euclidean statics general contain references to both 
single forces and to force couples. However, as is easy to verify by a simple 
calculation, the resultant of two such forces in homogeneous form yields a 
resultant of the form cgq, c ^ 0. This corresponds to the ideal line: just as a 
translation is integrated in the projective model as a rotation around a ideal 
point, a force couple appears in the projective model as a force carried by the 
ideal line, whose intensity is proportional to the distance between the two 
original lines, and to the magnitude of the original force vectors. 
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